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ABSTRACT 


DETERMINATION  OF  THE  ULTIMATE  CAPACITY  OF  ELASTOMERIC 
BEARINGS  UNDER  AXIAL  LOADING 

Elastomeric  bearings  are  used  as  base  isolation  systems  in  order  to  reduce  the  response 
of  buildings  to  earthquake  ground  motions.  In  order  to  facilitate  the  use  of  this 
technology,  the  National  Institute  of  Standards  and  Technology  (NIST)  has  published 
guidelines  for  testing  base  isolation  systems.  NIST  seeks  to  improve  the  current 
guidelines  for  determining  the  capacity  of  bearings  under  axial  loading.  Very  few 
full-scale  bearings  have  been  tested  to  failure  because  the  capacity  of  a typical  bearing 
is  often  beyond  the  loading  capacity  of  existing  testing  facilities.  This  research  has 
sought  to  determine  if  tests  to  failure  under  axial  loading  on  scale  models  can  be  used 
to  predict  the  axial  failure  load  of  a full-scale  bearing.  Tests  to  failure  were  conducted 
on  full-scale,  1/2-scale,  and  1/4-scale  bearings.  However,  there  were  dissimilarities  in 
steel  shim  thickness,  cover  layer,  and  the  central  alignment  hole.  The  stress-strain 
responses  were  similar,  but  not  identical.  Nonlinear  finite  element  analysis 
conclusively  demonstrates  that  similarity  in  shim  and  rubber  thickness  is  most 
important.  The  dissimilarities  in  the  cover  layer  and  central  alignment  holes, 
separately  considered,  had  less  effect  on  the  bearing  response  than  shim  thickness 
dissimilarities,  and  offset  each  other  in  this  particular  case.  Tests  indicate  that  the 
mechanical  properties  of  the  22-gage  shims  (0.762  mm  thick)  for  the  1/4-scale  bearing 
may  be  slightly  modified,  possibly  due  to  sandblasting  and/or  the  cold  rolling 
processing  of  the  steel. 

Key  Words:  building  technology;  compression  testing;  earthquake  engineering; 
elastomeric  bearings;  seismic  base  isolation;  ultimate  load  tests;  nonlinear  finite 
elements;  Valanis-Landel  Function;  rubber  modeling. 
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PREFACE 


Any  mention  of  commercial  items  serves  only  to  adequately  specify  the  analysis  and 
experimental  procedures.  In  no  case  does  such  identification  imply  recommendation 
or  endorsement  by  the  National  Institute  of  Standards  and  Technology,  nor  does  it 
imply  that  this  product  is  necessarily  the  best  available  for  the  purpose. 
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1.  INTRODUCTION 


Seismic  base  isolation  is  one  of  the  most  promising  concepts  in  seismic 
resistant  construction  to  come  of  age  in  this  century.  Most  of  the  energy  of  an 
earthquake  occurs  at  a period  very  near  the  fundamental  period  of  low  to  mid-rise 
buildings.  Base  isolation  reduces  the  response  of  the  structures  of  this  size  to 
earthquake  ground  motions  by  shifting  the  structure’s  fundamental  period  away  from 
the  dominant  period  of  the  ground  shaking.  In  an  effort  to  facilitate  the  use  of  this 
technology,  the  National  Institute  of  Standards  and  Technology  (NIST)  has  published 
Guidelines  for  Pre-Qualification,  Prototype  and  Quality  Control  Testing  of  Seismic 
Isolation  Systems  (NISTIR  5800,  1996). 

Base  isolation  systems  provide  the  lateral  flexibility  and  damping  necessary  for 
effective  isolation,  and  the  high  initial  stiffness  required  to  resist  wind  loads.  Isolation 
systems  are  broadly  grouped  into  three  categories:  elastomeric,  sliding,  and  hybrid.  A 
common  type  of  elastomeric  isolation  system  is  the  ordinary  laminated  bearing,  which 
consists  of  layers  of  elastomer  and  steel  bonded  under  high  temperature  and  pressure 
to  form  an  integral  bearing  that  is  free  of  joints  (elastomer  and  rubber  are  used 
interchangeably  throughout  this  report).  The  restoring  force  and  damping  mechanism 
are  combined  in  a single  device.  The  isolation  units  studied  in  this  research  are  of  the 
laminated  bearing  type.  Sliding  systems  rely  on  simple  Coulomb  friction  between  two 
surfaces,  and  hybrid  systems  can  incorporate  aspects  of  both  elastomeric  and  sliding 
types,  but  generally  use  independent  components  to  provide  the  restoring  force, 
damping,  wind  restraint  and  ultimate  restraint. 

Very  few  full-scale  isolation  units  have  been  tested  to  ultimate  failure  because 
the  bearing  capacity  of  a typical  bearing  is  often  beyond  the  loading  capacity  of 
existing  testing  facilities.  As  a result,  the  factor  of  safety  of  these  systems  is  unknown 
(NISTIR  5800,  Section  8.2,  1996).  A simple  way  around  this  problem  is  to  test  scale - 
model  versions  of  the  bearings.  If  the  material  and  geometric  properties  are  properly 
scaled,  then  similitude  relationships  can  be  used  to  predict  the  failure  load  of  a full- 
scale  bearing  based  on  tests  of  a scale  model  bearing.  However,  a perfectly 
geometrically  similar  scaled  model  is  difficult  to  make  because  some  dimensions  are 
predetermined  by  the  industry,  e.g.  standard  bearing  molds  are  frequently  not 
geometrically  similar,  resulting  in  dissimilar  diameters  and  central  cores.  Variation  of 
material  through  the  thickness,  if  any,  may  be  impossible  to  control.  The  major  goals 
of  this  study  are  to  establish  the  geometric  properties  that  must  be  similar  in  order  to 
build  scale  models  that  will  fail  at  the  same  stress  and  strain  as  a full-scale  bearing, 
and  to  determine  if  the  material  properties  of  the  steel  shims  change  significantly  in 
scale  models.  Four  laminated  bearings,  three  of  which  are  nearly  geometrically 
similar,  were  tested.  Finite  element  analyses  were  performed  to  investigate  the  effect 
of  geometric  dissimilarities  on  the  prediction  of  ultimate  load  from  scale  model  tests. 
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2.  UNIAXIAL  COMPRESSION  TESTING  ON  BEARINGS 

2.1  Description  of  Test  Specimens 

2.1.1  Dimensions 

Four  specimens  were  used  in  the  study.  The  specimens  are  designated  as  A 
and  B 1 through  B3.  Specimen  B1  is  a Full-Scale  (or  prototype)  elastomeric  bearing. 
Specimens  B2  and  B3  are  respectively  1/2-Scale  and  1/4-Scale  versions  of  Specimen 
B 1 . These  specimens  were  produced  by  the  same  manufacturer  with  the  goal  of 
minimizing  any  differences  in  material  properties  between  specimens.  Specimen  A is 
not  geometrically  similar  to  the  B-Bearings,  and  was  manufactured  by  a different 
company.  All  specimens  were  fabricated  using  materials  and  fabrication  procedures 
typical  of  those  used  in  the  manufacture  of  commercial  isolators.  Table  2.1  lists  the 
important  dimensions  of  each  bearing,  and  Figures  2.1  through  2.3.(b)  are 
dimensioned  drawings  of  Bearings  A,  B 1,  B2,  and  B3,  respectively. 

The  Full-Scale  Bearing  (Bl)  was  not  designed  for  a specific  application.  It 
was  sized  to  have  a slightly  lower  capacity  than  a typical  full-scale  isolator  in  order  to 
ensure  failure  in  the  testing  facilities  available  at  NIST. 

Even  though  strict  geometric  similarity  was  desired,  the  as-built  B-Bearings 
were  not  geometrically  similar.  The  dissimilarities  follow. 

( 1 ) Layer  Thicknesses  - For  geometric  similarity  to  hold,  each  shim  and 
rubber  layer  in  the  Full-Scale  Bearing  must  be  twice  as  thick  as  the  1/2-Scale  Bearing 
and  four  times  as  thick  as  1 /4-Scale  Bearing.  Due  to  a manufacturing  error  the  shims 
for  the  Full-Scale  Bearing  are  less  than  twice  the  thickness  of  the  shims  for  the  1/2- 
Scale  Bearing.  Geometric  similarity  in  shim  thickness  is  maintained  in  the  1/2-Scale 
and  1/4-Scale  Bearings,  and  in  rubber  thickness  for  all  B-Bearings. 

(2)  Bearing  Diameters  - The  outer  diameters  of  the  B-Bearings  could  not  be 
made  geometrically  similar  because  the  manufacturer  used  readily  available  molds. 
The  outer  diameters  of  the  bonded  areas,  which  correspond  to  the  shim  diameters,  are 
geometrically  similar.  Therefore,  the  cover  layers,  which  protect  the  shims  from 
environmental  exposure,  are  not  geometrically  similar. 

(3)  Center  Alignment  Hole  - The  center  alignment  holes  are  used  to  facilitate 
alignment  of  the  bearing  components  during  manufacture,  and  are  filled  with  rubber 
after  stacking  of  the  shim  and  rubber  layers.  They  are  not  geometrically  similar. 
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Table  2.1(a)  - Bearing  dimensions 


A 

Full-Scale 
B1  (mm) 

1/2-Scale 
B2  (mm) 

1 /4-Scale 
B3  (mm) 

Total  Height 

584 

304 

167 

109 

Total  Rubber  Layer  Thickness 

392 

191* 

95* 

48* 

Total  Shim  Thickness 

115.6 

37.2 

21.3* 

10.7* 

Thickness:  Each  End  Plate 

38.1 

38.1 

25.4 

25.4 

Table  2.1(b) 

- Bearing  dimensions 

A 

Full-Scale 

B1 

1/2-Scale 

B2 

1/4-Scale 

B3 

Rubber  Layers 

# of  Layers 

39 

15 

15 

15 

Thickness  (mm) 

10.1 

12.7* 

6.4* 

3.2* 

Shape  Factor** 

24.4 

15* 

15* 

15* 

Total  Diameter  (mm) 

1016 

800 

406 

229 

Bonded  Dia 

991 

762* 

381* 

191* 

(mm) 

Cover  Layer 
Thickness  (mm) 

2 x 12.5 

2x  19 

2 x 12.5 

2x  19 

* Geometrically  similar  dimension. 

**  Shape  factor  (S)  is  defined  as  the  bonded  plan  area  of  the  bearing  (one  side  only) 
divided  by  the  area  of  the  edge  of  a single  layer  of  elastomer  (S=d/4t),  where  d is  the 
bonded  diameter  and  t is  the  rubber  layer  thickness). 
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Table  2.1(b)  (continued)  - Bearing  dimensions 


A 

Full-Scale 

B1 

1/2-Scale 

B2 

1/4-Scale 

B3 

Steel  Shims 

No.  of  shims 

38 

14 

14 

14 

Thickness  (gage) 

11  gage* 

12  gage 

16  gage* 

22  gage* 

(mm) 

3.048* 

2.667 

1.524* 

0.762* 

Diameter  (mm) 

991 

762* 

381* 

191* 

Center  Alignment 

Hole  Diameter  t (mm) 

76.2 

50.1 

30.1 

30.1 

End  Plate 

Plan  Dimensions 
(mm) 

991  mm 
diameter 

876  mm  edge 
length 

508  mm  edge 
length 

305  mm  edge 
length 

* Geometrically  similar  dimension. 

t Holes  are  placed  in  shim  plates  to  facilitate  alignment  of  the  bearing  components 
during  manufacture. 

2.1.2  Rubber  Properties  Reported  by  Bearing  Manufacturers 

The  following  sections  present  the  results  of  ASTM  tests  on  the  rubber 
compounds,  as  reported  by  the  bearing  manufacturers.  Both  bearing  types  use  a high- 
damping  natural  rubber  vulcanizate. 

2.1.2.1  Rubber  for  B-Bearings 

All  natural  rubber  used  for  these  bearings  came  from  the  same  shipment.  The 
natural  rubber  is  mixed  with  various  chemicals,  but  due  to  the  size  of  the  rubber 
mixing  equipment,  three  batches  were  required.  Four  batches  were  mixed,  but  plant 
quality  control  rejected  the  first  batch  (the  rejected  batch  had  been  calendered  (rolled) 
to  the  wrong  thickness).  Each  batch  was  made  from  the  same  proportions  of  identical 
ingredients.  The  rubber  was  vulcanized  only  after  placing  the  shims  and  the  rubber  in 
alternating  layers  into  a mold. 

Two  bearings  of  each  size  were  manufactured;  one  bearing  of  each  size  was 
tested  in  this  program.  The  second  bearing  will  be  tested  to  failure  in  compression 
with  an  applied  lateral  load  in  a later  test  program.  Table  2.2  shows  the  batch  numbers 
for  each  bearing. 
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Table  2.2  - B-Bearings  - rubber  batches 


Bearing 

Ultimate  Compression 
No  Lateral  Load 

Ultimate  Compression  Lateral 
Load  Applied 

(Future  Test  Program) 

B1 

Batch  3 

Batch  2 

B2 

Batch  3 

Batch  3 

B3 

Batch  4 

Batch  4 

As  noted  above.  Batch  1 was  rejected  by  plant  quality  control,  not  because  its 
properties  differed  from  the  other  three  batches,  but  because  it  was  calendered  to  the 
wrong  thickness.  For  ASTM  testing,  four  of  the  six  tests  were  conducted  using 
specimens  made  from  a composite  batch,  which  consisted  of  equal  quantities  from 
each  of  the  four  batches.  The  results  of  ASTM  testing  are  shown  in  Table  2.3. 

2.1.2.2  Rubber  for  A-Bearing 

Test  results  of  the  rubber  used  in  the  particular  bearing  tested  were  not 
available.  The  manufacturer,  however,  reported  test  results  for  elastomers  with  the 
same  chemical  composition.  These  results  are  shown  in  Table  2.3. 

Table  2.3  - ASTM  test  results  for  elastomers 


Shear  Modulus 

(ASTMD 
4014-89, 
Annex  Al) 
(psi) 

Tensile 

Strength 

(ASTMD  412- 
92,  Method  A) 

(psi) 

Ultimate 

Elongation 

(ASTMD  412- 
92,  Method  A) 

(%) 

Compression 

Set 

(ASTMD  395- 
89,  Method  B) 

(%) 

B-Bearings 

Composite 

Batch 

128 

2406 

621 

20.75 

Batch  1 

N/A 

2425 

597 

N/A 

Batch  2 

N/A 

2537 

608 

N/A 

Batch  3 

N/A 

2491 

613 

N/A 

Batch  4 

N/A 

2589 

664 

N/A 

Bearing  A 
elastomer 

173 

2827-2891 

529  -579 

N/A 

N/A  denotes  that  tests  were  not  performed. 
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Table  2.3  (continued)  - ASTM  test  results  for  elastomers 


Shore  A Hardness 

(ASTM  D 2240-91) 

Peel  Strength 

(ASTM  D 429-81, 
Method  B) 

Ozone  Resistance 

(ASTM  D 1149-91)* 

B-Bearings 

Composite 

Batch 

65 

89  lb/in. 

100%  rubber  tear 

No  cracks  at 
magnification  7X 

Batch  1 

64 

N/A 

N/A 

Batch  2 

65 

N/A 

N/A 

Batch  3 

65 

N/A 

N/A 

Batch  4 

65 

N/A 

N/A 

Bearing  A 
elastomer 

63-65 

81-93  lb/in 

N/A 

* The  ozone  concentration  was  50  parts  per  hundred  million  (50pphm).  The 
specimens  were  stretched  to  20%  strain,  and  placed  for  100  hours  at  a temperature  of 
100°  F. 

2.2  Test  Machine  and  Bearing  Setup 

The  Universal  Testing  Machine  is  located  in  the  Large  Scale  Structural  Testing 
Facility  (Building  202)  at  the  National  Institute  of  Standards  and  Technology.  It  is 
managed  by  the  Building  and  Fire  Research  Laboratory.  The  test  machine  has  a 
maximum  compression  capacity  of  53.4  MN  (12  million  pounds).  See  Figure  2.4. 

The  machine  has  four  large  screw  columns  that  rise  to  approximately  23.8  m above 
the  ground.  The  load-head  is  connected  to  the  four  columns,  and  is  positioned  by 
moving  it  up  or  down  on  the  screw  columns.  Once  in  position,  the  head  is  locked  in 
place  so  that  the  head  does  not  move  relative  to  the  screw  columns  throughout  the  test. 
The  screw  columns  extend  to  about  7.0  m below  ground,  where  a hydraulic  cylinder 
assembly  is  located.  The  cylinder  is  approximately  2.1  m in  diameter  and  is  rigidly 
connected  to  the  four  screw  columns.  As  fluid  is  pumped  into  the  cylinder,  the 
cylinder  is  pulled  down,  which  in  turn  pulls  the  four  columns  down  simultaneously. 
Consequently,  the  load-head  moves  down,  pushing  the  head  against  the  specimen. 

The  maximum  stroke  of  the  machine  is  1 .5  m. 

For  the  A-Bearing  tests,  the  bearing  was  placed  on  two  large  cylindrical  steel 
plates.  These  were  in  turn  placed  on  a reinforced  concrete  platform  designed  to 
withstand  greater  than  53.4  MN.  See  Figure  2.5.  For  the  B-Bearing  tests,  the 
reinforced  concrete  platform  was  not  used.  Instead,  a third  cylindrical  piece  of  steel 
was  placed  on  top  of  the  two  steel  plates  used  in  the  A-Bearing  tests.  See  Figure  2.6. 
Also,  for  each  test  a smooth  steel  plate,  6.35  mm  (1/4  in)  thick,  was  placed  directly 
above  and  below  the  bearing  to  ensure  a smooth  contact  surface  and  to  protect  the  test 
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machine  load-head  from  damage  that  could  result  from  the  irregular  surface  on  the  top 
of  the  bearings. 

2.3  Control  of  Test  Machine 

The  test  machine  is  controlled  by  an  open-loop  control  system.  The  amount  of 
hydraulic  fluid  in  the  hydraulic  cylinder  assembly  determines  the  head  position,  and 
the  rate  of  flow  of  hydraulic  fluid  determines  the  head  speed.  To  control  the  head 
speed,  the  operator  manually  adjusts  the  rate  of  flow  of  hydraulic  fluid.  The  rate  of 
loading  is  indicated  by  the  speed  of  the  load  indicator  and  tells  the  operator  if  he  is 
allowing  too  much  or  too  little  fluid  to  move  the  head.  To  maintain  a constant  head 
velocity  (i.e.,  a constant  rate  of  flow  of  hydraulic  fluid),  the  operator  must  continually 
increase  the  valve  opening  to  compensate  for  a reduced  flow  rate  as  the  fluid  pressure 
rises.  However,  even  with  these  manual  adjustments,  head  speed  decreases  as  the  load 
approaches  the  machine’s  capacity  of  53.4  MN  (12  million  pounds). 

2.4  Test  Procedure 

As  described  in  NISTIR  5800,  Section  5.4.1,  each  bearing  was  subjected  to 
simple  axial  compression  under  zero  lateral  load.  Two  tests  were  conducted  on 
Bearing  A in  order  to  gain  experience  with  test  machine  control  and  data  acquisition 
prior  to  conducting  the  tests  on  the  geometrically  similar  bearings  Bl,  B2,  and  B3. 

It  is  the  usual  practice  in  testing  elastomeric  bearings  to  mechanically 
condition  the  elastomer,  since  the  initial  response  of  an  elastomeric  specimen  is  stiffer 
than  the  response  after  even  one  cycle  of  loading.  The  response  of  the  elastomer 
typically  stabilizes  after  a number  of  cycles;  however,  there  is  no  firm  guidance  on  the 
number  of  conditioning  cycles  needed  to  attain  equilibrium.  Researchers  of  bridge 
elastomeric  bearings  have  used  four  (4)  mechanical  conditioning  cycles  (Minor  and 
Egen,  1970).  For  dynamic  testing,  at  least  six  cycles  have  been  recommended  by 
Brown  (1986).  In  this  study,  it  was  decided  to  use  at  least  six  conditioning  cycles.  In 
the  first  test  on  Bearing  A,  seven  conditioning  cycles  to  13.79  MPa  (2000  psi) 
compressive  stress  were  applied.  In  the  second  test  on  the  A-Bearing  and  the 
subsequent  tests  on  the  B-Bearings,  six  conditioning  cycles  were  applied  to  13.79  MPa 
(2000  psi)  compressive  stress.  Six  conditioning  cycles  were  adequate  to  reach  an 
equilibrium  stress-strain  curve. 

2.5  Data  Acquistion 

Four  Linear  Variable  Displacement  Transducers  (LVDTs)  were  placed  evenly 
around  each  bearing  in  order  to  measure  vertical  displacement.  These  figures  were 
averaged  to  arrive  at  the  vertical  displacement.  The  LVDTs  were  not  attached  to  the 
bearings,  but  contacted  the  load  head  directly,  as  shown  in  Figure  2.6.  Therefore,  any 
misalignment  that  might  exist  between  the  end  plates  of  the  bearings  is  not 
measurable.  Displacement  readings  were  electronically  recorded  every  1/1 0th  of  a 
second. 
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2.6  Observations  and  Results  of  A and  B-Bearing  Tests 

2.6.1  Bearing  A - First  Test 

The  initial  test  on  the  A-Bearing  consisted  of  seven  mechanical  conditioning 
cycles  followed  by  five  load  cycles  to  the  machine  capacity  of  53.4  MN.  The  test  was 
conducted  at  the  maximum  head  speed,  which  was  measured  at  0.40  mm/sec. 
However,  Figure  2.7,  which  shows  a plot  of  head  speed  versus  load  for  a typical  load 
cycle,  clearly  shows  that  maintaining  a constant  head  velocity  was  not  possible.  The 
average  speed  over  this  cycle  was  0.30  mm/sec,  with  a standard  deviation  of  0.09 
mm/sec. 

During  the  mechanical  conditioning  cycles  to  13.79  MPa  (2000  psi),  no 
bulging  of  the  sides  of  the  bearing  was  observed.  During  the  load  cycles  to  the 
machine  maximum  of  53.4  MN,  slight  bulging  was  visible.  The  middle  third  of  the 
bearing,  as  viewed  from  the  West,  skewed  slightly  to  the  North.  There  was  a slight 
residual  bulge  and  skew  at  the  end  of  the  test. 

The  load  vs  displacement  plot  (Figure  2.8)  shows  a change  in  stiffness  at 
approximately  50  MN  compressive  force  (65  MPa)  and  a compressive  strain  of  0.042 
m/m  in  the  first  loading  cycle.  After  removal  of  load,  the  bearing  did  not  return  to  its 
pre-loaded  height,  thus  confirming  that  inelastic  deformation  took  place  during  the 
first  cycle.  It  is  reasonable  to  conclude  that  the  change  in  slope  in  the  first  loading 
cycle  was  due  to  yielding  of  the  steel  shims. 

Prior  to  initial  yielding,  the  curve  is  nearly  linear,  with  a tangent  stiffness  of 
1700  MPa.  After  the  initial  yielding,  the  tangent  stiffness  dropped  by  47%  to  900 
MPa.  The  curve  shows  a concave-upward  shape  for  the  second  through  fifth  load 
cycles.  The  stress-strain  response  of  pure  rubber  has  a concave-upward  shape  at 
higher  compressive  stresses,  and  this  supports  the  conclusion  that  the  steel  yielded  in 
the  initial  load  cycle  to  53.4  MN  - the  bearing’s  behavior  is  dominated  more  by  the 
rubber  than  the  steel  in  the  second  through  fifth  cycles. 

During  the  unloading  portion  of  the  first  cycle,  the  operator  inadvertently 
reversed  the  flow  of  hydraulic  fluid,  causing  a slight  increase  in  load  for  a moment. 
This  is  noticeable  in  the  plot. 

2.6.2  Bearing  A - Second  Test 

Rubber’s  stress-strain  behavior  is  a function  of  strain  rate  (Treloar,  1975).  The 
second  test  on  the  A-Bearing  was  conducted  primarily  to  investigate  the  influence  of 
strain  rate  on  bearing  response.  Four  sets  of  three  cycles  of  loading  were  applied,  each 
to  the  machine  load  capacity.  The  first  set  was  conducted  at  the  maximum  head 
speed,  and  the  second,  third  and  fourth  sets  at  50%  , 25%  , and  10%  of  the  maximum 
head  speed,  respectively  (the  maximum  head  speed  measured  in  these  tests  was  0.40 
rnm/s).  As  shown  in  Figure  2.9,  the  stress-strain  curves  for  each  set  are  virtually 
identical. 

The  second  test  was  performed  seven  days  after  the  first.  The  slight  bulge  and 
skew  of  the  bearing  in  the  northerly  direction  resulting  from  the  first  test  was  still 
visible  prior  to  the  second  test.  As  in  the  first  test,  the  bulging  and  lateral  deformation 
increased  when  the  bearing  was  loaded. 
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The  general  shape  of  the  curves  (Figure  2.9)  is  similar  to  the  curves  obtained 
in  the  first  test  after  yielding.  The  vertical  compressive  strain  in  the  first  load  cycle  of 
the  second  test  was  approximately  0.04  m/m  at  53.4  MN,  which  is  about  the  same  as 
the  vertical  strain  in  the  second  load  cycle  of  the  first  test  (the  first  cycle  after 
yielding),  with  the  permanent  set  due  to  yielding  removed.  The  plot  also  shows  a 
slight  increase  in  the  maximum  displacement  in  each  succeeding  load  cycle,  as  in  the 
first  test. 

The  plot  also  indicates  that  the  bearing  response  was  virtually  identical  at  head 
speeds  ranging  from  the  maximum  head  speed  down  to  approximately  10%  of  the 
head  speed.  Therefore,  for  the  limited  range  of  loading  rates  examined  in  this  study, 
bearing  response  to  compressive  loading  can  be  considered  independent  of  loading 
rate. 

2.6.3  Full-Scale  Bearing  (Bl)  Test 

For  each  of  the  B-Bearings  tests,  the  bearings  failed  as  a result  of  shim  rupture. 
Accordingly,  the  loading  history  consisted  of  six  mechanical  conditioning  cycles 
followed  by  loading  to  failure.  In  all  cases,  load  was  applied  well  beyond  initial  shim 
failure. 

2.6.3.1  Visual  Observations  and  Load  vs  Displacement  Curve  for  Full-Scale 
Bearing  (Bl)  Test 

Bulging  was  first  noticeable  at  approximately  65  MPa  compressive  stress  and 
12%  compressive  strain,  and  increased  sharply  between  80  and  90  MPa.  Bulging  was 
greatest  at  the  center,  and  varied  in  a parabolic  shape  from  the  top  to  the  bottom  of  the 
bearing.  Bulging  appeared  to  be  fairly  uniform  around  the  circumference  of  the 
bearing.  Failure  occurred  at  90  MPa  compressive  stress  and  a compressive  strain  of 
0.267  m/m.  Failure  was  indicated  by  a loss  of  strength  as  displacement  continued. 
This  phenomenon  was  accompanied  by  an  audible  rupture,  followed  shortly  thereafter 
by  a rapid  series  of  audible  ruptures  as  additional  shims  failed.  Each  audible  rupture 
was  followed  by  a sudden  drop  in  compressive  load.  Loading  was  stopped  when  the 
load  had  dropped  to  less  than  half  its  peak  value. 

The  stress-strain  curve  (Figure  2.10)  shows  initial  shim  yielding  at  a 
compressive  stress  of  about  42  MPa  and  a compressive  strain  of  about  0.06  m/m.  The 
shape  of  the  stress-strain  curve  up  to  this  point  is  concave-upward,  whereas  between 
initial  yield  and  the  first  shim  rupture  the  shape  is  concave-downward. 

After  the  initial  shim  rupture,  displacement  was  increased  at  the  same  rate.  As 
shims  failed,  the  load  dropped.  Small  load  recovery  occurred  between  shim  ruptures. 

2.6.3.2  Head  Velocity  and  Load  During  Failure  for  Full-Scale  Bearing  (Bl)  Test 

Figure  2.1 1 shows  load  and  head  velocity  versus  time  during  the  failure 
process.  Each  velocity  peak  corresponds  to  a shim  failure.  The  head  velocity  at 
failure  was  one  to  two  orders  of  magnitude  greater  than  the  head  velocity  prior  to 
failure.  This  is  consistent  with  the  visual  observation  that  the  head  movement  prior  to 
failure  was  imperceptible,  with  downward  jumps  of  the  load-head  observable  as  shims 
failed. 
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After  the  initial  shim  failure,  resistance  decreased  rapidly.  Five  shims  ruptured 
during  a period  of  nine  seconds,  and  the  resistance  decreased  from  40  MN  to  21  MN. 
The  bearing  then  recovered  to  22  MN,  whereupon  another  series  of  ruptures  followed. 
This  pattern  was  repeated  once  more  before  the  test  was  terminated. 

2.6.4  1/2-Scale  Bearing  (B2)  Test 

2.6.4.1  Visual  Observations  and  Load  vs  Displacement  Curve  for  1/2-Scale 
Bearing  (B2)  Test 

During  loading  to  failure,  bulging  became  noticeable  at  a compressive  stress  of 
about  40  MPa  and  a compressive  strain  of  about  0.052  m/m.  The  bulging  was 
asymmetric:  the  bearing  bulged  significantly  on  the  south  and  southwest  side,  and 
bulged  very  little  on  the  north  and  west  side.  This  bulging  pattern  continued  up 
through  failure.  Failure  occurred  at  88.7  MPa  compressive  stress  and  a compressive 
strain  of  0.227  m/m,  and  was  initiated  by  a shim  rupture  at  or  near  mid-height, 
followed  by  a series  of  ruptures  successively  further  away  from  mid-height. 

Figure  2.12  shows  initial  shim  yielding  at  a compressive  stress  of  about  45 
MPa  and  a compressive  strain  of  about  0.06  m/m.  Similar  to  the  Full-Scale  Bearing 
(B 1 ),  the  shape  of  the  stress-strain  curve  up  to  this  point  is  concave-upward,  whereas 
between  initial  yield  and  the  first  shim  rupture  the  shape  is  concave-downward. 

2.6.4.2  Head  Velocity  and  Load  During  Failure  for  1/2-Scale  Bearing  (B2)  Test 

Figure  2.13  shows  load  and  head  velocity  versus  time  during  the  failure 
process.  Each  velocity  peak  corresponds  to  a shim  failure.  The  head  velocity  at 
failure  was  one  to  two  orders  of  magnitude  greater  than  the  velocity  prior  to  failure. 
This  is  consistent  with  the  visual  observation  that  the  head  movement  prior  to  failure 
was  imperceptible,  with  downward  jumps  of  the  load-head  observable  as  shims  failed. 

After  initial  failure,  load  fell  quickly.  Over  fifteen  seconds,  there  were  at  least 
five  shim  ruptures,  and  the  resistance  decreased  from  10.7  MN  to  5.8  MN.  The 
bearing  then  recovered  to  6.8  MN,  whereupon  there  were  another  series  of  ruptures. 

2.6.5  1/4-Scale  Bearing  (B3)  Test 

2.6.5.1  Visual  Observations  and  Load  vs  Displacement  Curve  for  1/4-Scale 
Bearing  (B3)  Test 

Bulging  was  first  noticeable  at  a compressive  stress  of  about  20  MPa  and  a 
compressive  strain  of  about  0.067  m/m.  Bulging  increased  steadily  to  failure,  at  a 
compressive  stress  96.7  MPa  and  a compressive  strain  of  0.281  m/m.  Bulging 
appeared  to  be  uniform  around  the  circumference  of  the  bearing.  Bulging  was  greatest 
at  the  mid-height  and  the  shape  was  parabolic  when  viewed  from  the  side.  Initial  shim 
failure  was  not  audible.  It  appeared  that  shim  rupture  was  at  or  near  the  mid-height  of 
the  bearing,  with  subsequent  failures  occurring  successively  further  from  mid-height. 

In  this  test,  load  application  was  continued  considerably  beyond  failure  of  all 
shims.  The  load  carrying  capacity  eventually  increased  because  after  the  failure  of  all 
shims,  the  bearing  essentially  acted  as  a single,  unbonded  layer  of  elastomer  (even 
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bond  to  the  end  plates  was  eventually  lost).  By  the  end  of  the  test  compressive  strains 
exceeded  0.75  m/m. 

Figure  2.14  shows  initial  shim  yielding  at  a compressive  stress  of  about  52 
MPa  and  a compressive  strain  of  about  0. 1 1 m/m.  Similar  to  the  Full  and  1/2-Scale 
Bearings,  the  shape  of  the  stress-strain  curve  up  to  this  point  is  concave-upward, 
whereas  between  initial  yield  and  the  first  shim  rupture  the  shape  is  concave - 
downward. 

2.6.5.2  Head  Velocity  and  Load  During  Failure  for  1/4-Scale  Bearing  (B3)  Test 

Figure  2.15  shows  load  and  head  velocity  versus  time  during  the  failure 
process.  Each  velocity  peak  corresponds  to  a shim  failure.  The  head  velocity  at  failure 
was  one  order  of  magnitude  greater  than  the  velocity  prior  to  failure.  This  is 
consistent  with  the  visual  observation  that  the  head  movement  prior  to  failure  was 
imperceptible,  with  downward  jumps  of  the  load-head  observable  as  shims  failed. 

After  initial  failure,  load  fell  quickly.  Over  thirty  seconds,  there  were  at  least 
five  shim  ruptures,  and  the  resistance  decreased  from  2.75  MN  to  1.70  MN.  The 
bearing  then  recovered  to  1 .85  MN,  whereupon  there  were  another  series  of  ruptures. 

2.6.6  Comparisons  among  A and  B-Bearing  Tests 

Figure  2.16  shows  the  stress-strain  curves  of  all  specimens  on  the  same  plot 
(the  first  cycle  of  loading  to  53.4  MN  for  the  A-Bearing  test  is  shown).  Note  the 
similarity  between  the  three  curves  for  the  three  B-Bearings,  in  spite  of  the  fact  that 
the  bearings  were  not  strictly  geometrically  similar  (c.f.  Table  2.1).  All  curves  of  the 
B-Bearings  have  similar  slopes,  except  at  low  stress  values.  Differences  in  the  initial 
stiffnesses  may  be  due  in  part  to  misalignment  of  the  load-head  and  endplate. 
Measurement  of  any  misalignment  is  not  possible,  since,  as  shown  in  Figure  2.6,  and 
as  discussed  in  Section  2.5,  the  LVDTs  were  not  attached  to  the  bearings,  but  instead 
contacted  the  load  head  directly. 

The  approximate  tangent  stiffness  of  the  A-Bearing  is  1700  MPa  just  prior  to 
yield.  The  corresponding  stiffness  for  the  B-Bearings  is  930  MPa.  Bearing  A had  a 
rubber  layer  shape  factor  of  24.4,  whereas  the  B-Bearings  had  a shape  factor  of  15.  A 
higher  shape  factor  results  in  a higher  compressive  stiffness  for  the  A-Bearing. 

Similar  observations  have  been  reported  by  researchers  of  bridge  elastomeric  bearings 
(Minor  and  Egen,  1970;  Stanton  and  Roeder,  1982;  Roeder  etal.,  1987). 

If  the  lower  initial  stiffnesses  of  the  Full-Scale  and  1 /4-Scale  Bearings  were 
indeed  due  to  plate  and  head  misalignment,  the  stress-strain  curves  can  be  shifted 
horizontally.  In  Figure  2.17,  the  1/4-Scale  stress-strain  and  Full-Scale  stress-strain 
curves  were  shifted  to  the  left  by  strains  of  .037  and  .005  respectively,  in  order  to 
clearly  show  the  nearly  identical  slope  of  the  stress-strain  curves  just  prior  to  initial 
shim  yielding.  However,  the  stress-strain  curves,  while  similar  in  shape,  are  clearly 
not  identical.  For  example,  the  maximum  stress  at  failure  for  the  Full-Scale  Bearing  is 
90  MPa;  for  the  1/2-Scale  Bearing  the  maximum  stress  is  88.7  MPa,  or  1.4%  smaller, 
and  for  the  1/4-Scale  Bearing  the  maximum  stress  is  96.7  MPa,  or  7.4%  larger.  Also, 
the  strain  at  failure  for  the  Full-Scale  Bearing  is  0.267,  but  is 
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0.227  (-13.0%  ) for  the  1/2-Scale  Bearing,  and  0.281  (+6.5%)  for  the  1/4-Scale 
Bearing.  Finally,  there  is  a clear  trend  toward  higher  stress  and  strain  at  the  point  of 
initial  yielding  as  we  move  from  the  Full-Scale  Bearing  down  to  the  1/4-Scale 
Bearing. 

2.7  Measurement  Accuracy  and  Variability 

The  accuracy  of  the  Universal  Testing  Machine  machine  is  established  through 
a calibration  using  known  dead  loads,  and  was  calibrated  according  to  ASTM  E4, 
which  requires  that  the  machine  reading  be  within  ± 1 % of  the  true  reading.  The 
machine  has  met  these  requirements.  This  is  equivalent  to  a relative  expanded 
uncertainty  (k=2)  of  1%,  and  is  small  compared  to  specimen  variability. 

The  bearing  testing  program  consisted  of  only  one  test  of  each  size  of  bearing, 
and  the  stresses  at  bearing  failure,  from  lowest  to  highest,  are  88.7,  90.0,  and  96.7 
MPa  for  the  1/2-Scale,  Full-Scale,  and  1/4-Scale  Bearings,  respectively.  The  failure 
stresses  range  from  -1.4%  to  +7.4%  of  the  median  value,  so  the  variability  in  these 
measurements  can  be  attributed  to  factors  outside  the  resolution  accuracy  of  the  test 
machine. 

Sources  of  variation  of  failure  stresses  include  (1)  deviations  from  perfect 
geometric  similarity,  (2)  differences  in  the  steel  and  rubber  properties,  and  (3) 
imperfect  manufacture.  In  Section  6.5,  the  effect  of  geometric  dissimilarities  on 
bearing  response,  assuming  identical  steel  and  rubber  properties,  is  demonstrated 
through  finite  element  analysis.  Results  of  steel  tensile  coupon  tests,  summarized  in 
Section  5.3,  show  higher  values  for  the  modulus  of  elasticity  and  stresses  in  the  plastic 
range  for  the  1/4-Scale  shims  (22  gage),  and  the  effect  on  bearing  response  is 
discussed  in  Section  6.5.  The  effects  of  rubber  property  variability  have  not  been 
determined  in  this  study  (c.f.  Section  3.3.5),  nor  have  the  effects  of  imperfect 
manufacture. 
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Figure  2.1  - Bearing  A dimensions 
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Figure  2.2  - Full-Scale  Bearing  (Bl)  dimensions 
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(a)  1 /2-Scale  Bearing  (B2)  dimensions 


(b)  1/4-Scale  Bearing  (B3)  dimensions 

Figure  2.3  - Bearing  dimensions  for  1/2-Scale  (B2)  and  1/4-Scale  (B3)  Bearings 
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Figure  2.4  - Test  machine 
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Figure  2.5  - Setup  for  Bearing  A tests 
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Figure  2.6  - Setup  for  geometrically  similar  bearings  (B)  tests 
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Figure  2.7  - Load-head  speed  vs  load  for  Bearing  A,  initial  test  (2nd  cycle) 
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Figure  2.8  - Load  vs  displacement  plot  for  Bearing  A,  first  test 
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Figure  2.9  - Load  vs  displacement  plot  for  Bearing  A,  second 
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Figure  2.10  - Load  vs  displacement  plot  for  Full-Scale  Bearing  (Bl) 
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Figure  2.1 1 - Head  velocity  and  load  during  shim  failures  for  Full-Scale  Bearing  (Bl) 
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Figure  2.12  - Load  vs  displacement  plot  for  1/2-ScaleBearing  (B2) 
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Figure  2.13  - Head  velocity  and  load  during  shim  failures  for  1/2-Scale  Bearing  (B2) 
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Figure  2.14  - Load  vs  displacement  plot  for  1/4-Scale  Bearing  (B3) 
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Figure  2.15  - Head  velocity  and  load  during  shim  failures  for  1/4-Scale  Bearing  (B3) 
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Figure  2.16  - Comparison  of  stress-strain  curves  for  all  bearings 
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Figure  2.17  - Comparison  of  adjusted  stress-strain  curves  for  B-Bearings 
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3.  DETERMINATION  OF  RUBBER  STRESS-STRAIN  PROPERTIES  FOR 
FINITE  ELEMENT  MODELING 

3.1  Background  for  Finite  Element  Modeling 

Finite  element  analyses  have  been  conducted  on  the  B-Bearings  to  investigate 
the  effect  of  the  geometric  dissimilarities  on  the  performance  of  the  scale-model 
bearings,  and  their  abilities  to  predict  the  ultimate  load  capacity. 

As  noted  in  Section  2.3.3,  each  bearing  was  subjected  to  at  least  six 
mechanical  conditioning  cycles  to  13.79  MPa  before  they  were  loaded  to  failure. 
During  these  cycles,  the  bearing  deformation  was  not  elastic.  Figure  3.1  shows  a 
tensile  load-displacement  plot  of  a dumbbell  rubber  specimen,  where  there  were  three 
cycles  to  29.4  N.  The  plot  shows  a reduction  in  load  at  each  strain  level  for  the  second 
and  subsequent  cycles,  and  the  curve  takes  the  shape  of  the  initial  loading  cycle  when 
the  extension  is  continued  beyond  the  extension  reached  in  each  cycle.  A similar 
response  is  present  in  compressive  loading.  In  the  actual  bearing  tests,  the 
conditioning  cycles  took  the  bearings  to  only  10%  of  the  failure  strain  and  15%  of  the 
failure  stress.  Modeling  the  response  of  the  bearing  to  the  conditioning  cycles  would 
have  no  impact  on  capturing  the  initial  yielding  of  the  shims,  nor  on  the  response  of 
the  bearing  in  the  failure  range.  Therefore,  this  effect  was  not  modeled. 

The  response  of  rubber  is  also  a function  of  strain  rate  and  temperature 
(Treloar,  1975).  However,  we  are  not  considering  these  effects  in  this  analysis. 
Therefore,  a hyperelastic  model  was  employed  to  characterize  the  behavior  of  the 
rubber.  In  Section  3.2,  the  constitutive  theory  and  finite  element  modeling  of  rubber 
samples  are  discussed  in  detail.  In  Section  3.3,  the  experimental  testing  methods  and 
experimental  results  are  presented  (uniaxial  tension,  uniaxial  compression,  and 
volumetric  testing).  In  Section  3.4,  extension  of  the  experimental  data  to  other  modes 
of  deformation  is  presented  in  an  attempt  to  improve  the  characterization  of  the  strain- 
energy  function.  The  calibration  of  the  rubber  model,  e.g.  how  many  terms  are  to  be 
included  in  the  strain  energy  function,  and  whether  the  extension  of  the  data  to  other 
deformation  modes  improved  the  modeling,  is  discussed  in  Section  4. 

Steel  samples  of  12,  16,  and  22  gage  thicknesses  (c.f.  Table  2.1)  were  tested  in 
the  uniaxial  tensile  mode.  In  Section  5,  the  experimental  steel  testing  and  modeling 
will  be  discussed  in  detail.  Finally,  in  Section  6,  the  results  of  the  finite  element 
modeling  of  the  bearings  are  presented. 

The  finite  element  code  ABAQUS  was  used  to  model  the  B-Bearings 
responses  to  uniaxial  compressive  loading.  ABAQUS  is  a general  purpose  nonlinear 
finite  element  program  that  has  constitutive  models  which  are  appropriate  to  the 
analysis.  A hyperelastic  model  is  available  in  order  to  model  rubber-like  compounds, 
including  the  ability  to  model  slight  compressibility  of  rubber.  The  model  requires 
nominal  stress-strain  data  of  the  rubber.  For  modeling  the  steel,  an  elastic-plastic 
model  with  strain-hardening  was  needed,  since  uniaxial  tests  of  steel  samples 
indicated  this  type  of  response.  ABAQUS  provides  a model  with  these  characteristics, 
and  it  requires  true  (Cauchy)  stress  and  logarithmic  strain.  Therefore,  experimental 
steel  data  were  obtained  using  these  measures. 

The  input  file  for  ABAQUS  is  divided  into  two  parts  - (1)  the  model  data  and 
(2)  the  history  data.  The  model  data  define  the  geometry  of  the  model,  the  type  or 
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types  of  elements  used,  and  the  definitions  of  the  materials.  The  history  data  describes 
how  the  loads  and  displacements  are  prescribed.  ABAQUS  allows  multiple  steps  to 
be  defined,  where  each  step  may  be  a particular  set  of  loads  or  displacements,  and  the 
initial  condition  of  a step  is  the  solution  of  the  previous  step. 

To  model  the  simple  monotonic  loading  of  the  tests,  only  one  step  was 
required  in  the  history  data.  The  history  of  the  step  was  defined  by  specifying  the 
vertical  downward  displacement  of  the  top  nodes  to  be  30  mm  for  the  prototype,  15 
mm  for  the  1/2-Scale  Bearing,  and  7.5  mm  for  the  1/4-Scale  Bearing,  which 
corresponds  to  an  engineering  compressive  strain  of  0.314  m/m.  In  the  experimental 
tests,  the  compressive  strain  at  failure  ranged  from  0.227  m/m  for  the  1/2-Scale 
Bearing  to  0.281  m/m  for  the  1/4-Scale  Bearing.  A single  node  along  the  top  was 
designated  as  the  master  node,  and  the  vertical  displacement  of  all  other  nodes  along 
the  top  were  tied  to  the  vertical  displacement  of  the  master  node.  This  was 
convenient,  since  by  defining  the  displacements  this  way,  the  reaction  forces  for  all 
top  nodes  were  summed  at  the  master  node. 

Since  the  rate  of  loading  in  the  experiments  was  very  slow,  inertial  effects 
were  negligible.  Therefore,  a static  stress  analysis  was  performed.  This  type  of 
analysis  is  specified  in  ABAQUS  by  using  the  STATIC  keyword.  In  addition,  an  arc - 
length  control  method  was  used.  In  arc-length  control,  the  distance  traveled  along  the 
equilibrium  path  is  specified,  and  the  actual  load  value  is  unknown  at  the  outset. 
Without  the  use  of  the  arc-length  control  method,  convergence  could  not  be  obtained 
beyond  approximately  4%  of  the  desired  vertical  displacement.  Successful 
convergence  was  obtained  for  the  entire  vertical  displacement  with  the  use  of  arc- 
length  control.  The  arc-length  control  method  is  specified  in  ABAQUS  by  selecting 
the  keyword  RIKS.  For  further  information  on  these  topics,  consult  the 
ABAQUS/Standard  Users  Manual  (1995). 

3.2  Theoretical  Basis  for  Constitutive  Modeling  of  Rubber 

Rubber  consists  of  many  interwoven,  long  chain  polymers  that  initially  (prior 
to  loading)  are  randomly  oriented;  therefore,  rubber  can  be  assumed  to  be  isotropic. 
Also,  rubber  exhibits  a nonlinearly  elastic  response  up  to  large  strains,  i.e.  a loaded 
specimen  will  return  nearly  to  its  unloaded  state  after  being  stretched  or  compressed  a 
great  amount  (for  rubber,  the  tensile  strains  can  exceed  600%,  and  the  compressive 
strains  can  exceed  85%).  Consequently,  a hyperelastic  model  was  used  to  characterize 
the  rubber  behavior. 

A hyperelastic  model  associates  a state  of  strain  to  a unique  state  of  stress 
(“hyper”  means  that  the  model  is  elastic  over  a large  strain  range).  The  strain  energy 
potential  can  be  defined  at  each  point  in  the  material  based  on  this  stress-strain 
relationship,  and  stresses  can  be  found  directly  by  taking  appropriate  derivatives  of  the 
strain  energy  potential.  Since  the  material  is  assumed  to  be  isotropic,  the  strain  energy 

potential  can  be  written  as  a function  of  the  strain  invariants,  U = £/(/,, /2,/3),  where 

/3  is  the  square  of  the  volume  ratio  = j)  and  is  equal  to  1 for  a perfectly 

incompressible  material  (see  Appendix  A for  a description  of  the  invariants  of  the 
deformation  tensor  used  in  ABAQUS). 
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In  general,  the  total  stress  at  each  point  in  a body  can  be  found  if  the 
displacements  are  known.  However,  if  a material  is  perfectly  incompressible,  the 


hydrostatic  stress 


h 1 ^22  *33 


, where  t represents  true  stress,  cannot  be  found  from 


the  displacements,  since  the  application  of  a hydrostatic  stress  results  in  no  deformation. 
“Mixed”  formulations  have  successfully  dealt  with  this  problem  (Peng  and  Chang, 

1997;  Chen  and  Pan,  1996).  In  a mixed  formulation  for  a perfectly  incompressible 
material,  the  internal  energy  is  augmented  by  adding  the  term  p(J  - 1),  where  p is  a 
Lagrange  multiplier  (the  hydrostatic  stress)  introduced  to  impose  the  constraint 
7-1  = 0 (7-1  is  the  volumetric  strain).  This  allows  the  hydrostatic  stress  to  be 
approximated  directly,  independently  of  the  displacements.  The  stress  that  is  derived 
from  the  displacements  is  the  deviatoric  stress.  The  sum  of  the  hydrostatic  and 
deviatoric  stress  tensors  gives  the  total  stress  tensor. 


Good  results  have  been  obtained  when  modeling  solid  rubber  as  an 
incompressible  material.  A Poisson's  Ratio  of  0.5  is  equivalent  to  an  incompressible 
material.  For  solid  rubbers,  Poisson’s  Ratio  is  very  close  to  0.5,  usually  0.498  and 
above.  For  the  rubber  in  the  B-bearings,  it  is  estimated  to  range  from  0.4994  to  0.4999, 
so  it  is  very  close  to  an  incompressible  material.  This  range  of  values  was  found  from 
If  E\ 

the  equation  v = — ^1  - , which  requires  estimates  of  Young’s  modulus  E and  bulk 

modulus  K.  The  tests  to  determine  these  parameters  were  performed  on  rubber  samples 
taken  from  rubber  batches  used  to  produce  the  B-Bearings,  and  are  discussed  in  detail  in 
Section  3.3. 


The  lower  value  of  0.4994  was  found  as  follows:  Young’s  modulus  E was 
determined  by  averaging  the  moduli  from  uniaxial  compressive  and  tensile  tests,  and 
gave  a value  of  approximately  10.77  MPa.  The  bulk  modulus  K was  determined  from  a 
volumetric  compression  test  to  be  approximately  3250  MPa  (the  slope  of  the  Pressure 
versus  Volume  Ratio  curve,  Figure  3.13).  These  values  of  E and  K give  a Poisson’s 
ratio  of  0.4994. 


The  upper  value  of  0.4999  was  found  by  using  the  value  of  the  shear  modulus 
(G)  reported  by  the  bearing  manufacturer.  The  shear  modulus  was  determined  to  be 
128  psi  (0.8825  MPa).  The  bearing  manufacturer  used  the  double  shear  test  of  Annex 
A1  of  ASTM  D 4014-89  (c.f.  Table  2.3).  Young’s  modulus  E was  found  from  the 


where  a value  for  the  Poisson’s  ratio  of  0.4900  was  assumed. 


equation  G = —7—  . 

4 2(1 +v) 

The  resulting  value  of  Young’s  modulus,  and  the  value  of  the  bulk  modulus  discussed 

If  E \ 

above,  were  used  to  determine  Poisson’s  ratio  from  v = — 1 . This  value  was 

2 \ 3 K f 

compared  to  the  assumed  value  of  0.4900.  A second  iteration  showed  Poisson’s  ratio  to 
be  0.4999. 


Due  to  the  near  incompressibility  of  rubber  (v  =0.5),  when  loaded  it  will  tend 
to  maintain  its  volume  and  squeeze  into  any  free  space  available.  However,  there  isn’t 
much  free  space  available  for  movement  of  the  rubber  in  our  bearings  (Shape  Factor  = 
15),  except  at  the  exterior  of  the  bearing.  Since  the  rubber  is  highly  confined,  the  slight 
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compressibility  of  the  rubber  may  be  significant  in  obtaining  accurate  predictions  of  the 
stiffness  of  the  bearings,  and  therefore  could  affect  the  analytical  prediction  of  the 
deformation  and  ultimate  bearing  load  (Roeder  et  al.,  1987). 

As  noted  above,  the  total  stress  is  derivable  from  the  displacements,  assuming 
the  material  is  not  incompressible.  However,  the  near  incompressibility  of  solid  rubber 
means  that  the  stiffness  matrix  becomes  nearly  singular.  This  problem  is  compounded 
by  the  approximations  inherent  in  any  finite  element  approach,  for  example,  the 
calculation  of  the  strains  involves  taking  derivatives  of  the  displacement  numerically, 
resulting  in  more  error.  It  has  been  found  that  a mixed  formulation  can  be  utilized 
effectively  in  this  situation  (Peng  and  Chang,  1997;  Chen  and  Pan,  1996).  The  general 
approach  is  to  decouple  the  strain  energy  into  a purely  deviatoric  portion  and  a purely 
volumetric  portion.  The  deviatoric  portion  is  written  using  revised  invariants  that 
remove  any  effect  due  to  volume  change.  The  strain  energy  is  written  as 


u = f/(7„72,/3)  = (3.1) 

where  the  bars  over  the  first  two  strain  invariants  indicate  removal  of  volume  changes 
(see  Appendix  A). 

The  decoupling  of  the  deviatoric  and  volumetric  strain  energy  is  in  general  not 
valid,  but  it  is  required  for  the  use  of  the  mixed  approach.  It  has  been  shown  that  a 
constant  bulk  modulus  is  required  for  the  mixed  formulation  to  be  a valid  approach 
(Sussman  and  Bathe,  1987).  A constant  bulk  modulus  is  equivalent  to  linearity  between 
pressure  and  volume  ratio.  Figure  3.1 1 shows  a plot  of  pressure  versus  volume  ratio  for 
our  rubber,  as  determined  in  a volumetric  compression  test  The  curve  is  initially  flat, 
but  this  is  probably  due  the  rubber  specimen  not  fitting  exactly  into  the  testing 
apparatus,  thereby  allowing  the  rubber  to  flow  into  the  free  space.  However,  the  curve 
is  nearly  linear  beyond  that  region,  and  on  this  basis,  the  decoupling  of  the  strain  energy 
into  deviatoric  and  volumetric  portions,  as  shown  in  Equation  (3.1),  is  assumed  valid 
for  the  rubber  tested,  and  therefore  the  mixed  approach  is  an  appropriate  model. 

The  ABAQUS  hyperelastic  model  follows  the  mixed  approach,  and  was  used  in 
the  finite  element  analysis.  Two  forms  of  the  strain  energy  function  are  provided  in 
ABAQUS.  The  first  is  the  Polynomial  Strain  Energy  Function: 

u-  2c,/7,-3M-3)J  + |;-^(./,,-l)2i.  (3.2) 

i+j-l  i=l  L>i 

The  first  summation  is  the  contribution  due  to  deviatoric  effects,  and  the  second 
summation  is  the  contribution  due  to  volumetric  effects.  The  Cy  and  Z)  are  parameters 

that  are  found  from  the  rubber  test  data.  lx  and  J2  are  the  first  and  second  invariants  of 
the  Cauchy-Green  Deformation  Tensor  with  the  volume  change  removed,  and  Jel  is  the 
ratio  of  the  current  volume  to  the  original  volume  excluding  thermal  effects  (see 
Appendix  A).  For  N=  1,  the  deviatoric  contribution  to  the  strain  energy  is  the  Mooney- 
Rivlin  function,  which  is  U = C10  (7j  - 3)  + C01(J2  - 3). 

The  strain  energy  potential  function  of  an  isotropic  elastic  material  can  also  be 
written  in  terms  of  the  stretch  (extension)  ratios  A.  According  to  Valanis  and  Landel 


30 


(1967),  “A  necessary  and  sufficient  condition  that  U will  pertain  to  an  isotropic  elastic 
material  is  that  it  merely  be  a symmetric  function  of  the  extension  ratios  Xa . Thus,  an 
alternate  and  equally  general  form  of  U could  be  U = U(  ApA^Aj) Furthermore, 
Valanis  & Landel  postulate  that  the  strain  energy  function  for  rubber  can  be  written  as 
an  additive  and  separable  function  of  the  stretch  ratios  as  U = u(  A,)  + + m(A 3) 

(«(•)  is  termed  the  Valanis-Landel  function). 

The  second  form  for  the  strain  energy  potential  used  in  ABAQUS  is  the  Ogden 
Strain  Energy  Potential.  The  contribution  due  to  deviatoric  effects  is  written  as  a 
function  of  the  stretch  ratios,  and  is  functionally  identical  to  the  Valanis-Landel  form 
shown  above.  The  contribution  due  to  volumetric  effects  is  the  same  as  that  used  for 
the  Polynomial  Strain  Energy  Function: 


(3.3) 


In  Equation  (3.3),  the  ju.  and  a,  are  parameters  that  are  found  from  the  rubber 


test  data.  The  principal  stretch  ratios  \ have  the  volume  change  removed  (see 


Appendix  A).  To  determine  the  deviatoric  parameters  C ip  and  a„  experimental  data 
can  be  fitted  to  one  or  more  of  three  deformation  modes;  uniaxial  (tension  or 
compression),  equibiaxial  (tension  or  compression),  and  planar  (tension  or 
compression)  (ABAQUS  1995a).  For  a right-handed  Cartesian  coordinate  system  with 
axes  xl,  x2,  and  x3,  the  uniaxial  deformation  mode  corresponds  to  an  application  of 
force  along,  e.g.,  the  longitudinal  axis  of  the  specimen.  For  the  purpose  of  explanation 
here,  we  denote  this  axis  xl.  Then  the  specimen  is  free  to  expand  or  contract  along  axes 
x2  and  x3.  The  equibiaxial  deformation  mode  corresponds  to  equal  stretches  along  the 
xl  and  x2  axes,  with  unrestrained  movement  along  the  x3  axis,  and  planar  deformation 
corresponds  to  a stretch  along  the  xl  axis,  with  no  stretch  allowed  along  the  x2  axis, 
and  unrestrained  movement  along  the  x3  axis.  Each  of  these  modes  is  an  example  of  a 
deformation  in  which  the  directions  of  principal  strain  do  not  change.  Each  mode  can 
be  written  in  terms  of  a single  stretch,  which  allows  the  stress-strain  relationship  to  be 
easily  measured. 

For  an  incompressible  material,  uniaxial  compression  is  equivalent  to 
equibiaxial  tension,  uniaxial  tension  is  equivalent  to  equibiaxial  compression,  and 
planar  tension  is  equivalent  to  planar  compression  (ABAQUS  1995a).  For  this  analysis, 
experimental  testing  was  performed  in  uniaxial  tension  and  uniaxial  compression,  and  is 
described  in  detail  in  Sections  3.3.1  and  3.3.2.  It  was  felt  necessary  to  include  both 
tension  and  compression  data.  The  exclusion  of  compression  data  may  result  in  a strain 
energy  function  that  does  not  capture  much  of  the  strain  energy  dependence  on  /2.  A 
planar  test  would  have  been  done  had  we  had  the  appropriate  equipment,  and  may  have 
captured  dependence  on  the  strain  invariants  that  was  missed  by  the  uniaxial  tests.  In 
lieu  of  an  actual  planar  test,  we  calculated  the  stress-strain  response  in  planar  tension 
using  the  Valanis-Landel  function  (c.f.  Section  3.4).  In  order  to  determine  the 
volumetric  parameter  D, , which  is  needed  to  model  compressibility,  an  experimental 
volumetric  test  was  performed.  The  volumetric  test  and  results  are  described  in  detail  in 
Section  3.3.3. 
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3.3  Experiments  Performed  on  Rubber  for  Finite  Element  Modeling 
3.3.1  Rubber  Uniaxial  Tension  Testing 

The  uniaxial  tension  tests  were  carried  to  engineering  tensile  strains  of  6.00 
m/m,  very  near  the  ultimate  elongation  of  the  rubber,  in  order  to  ensure  that  the  strain 
energy  function  would  adequately  represent  the  severe  engineering  tensile  stresses  and 
strains  expected  in  the  finite  element  analysis. 

The  bearing  manufacturer  supplied  two  vulcanized  sheets  of  rubber  for  each  of 
the  four  batches  and  the  composite  batch,  which  consisted  of  equal  proportions  of  the 
four  batches.  As  discussed  in  Section  2. 1.2.1,  the  bearings  were  manufactured  from 
three  batches  of  rubber.  A fourth  batch  was  rejected  because  it  was  calendered  (rolled) 
to  the  wrong  thickness.  Therefore,  a sample  from  the  composite  batch  is  representative 
of  each  of  the  batches.  Each  sheet  was  prepared  in  accordance  with  ASTM  D 3 182-89. 

The  calendering  of  the  rubber  sheets  may  have  a measurable  influence  on  tensile 
properties  since  anisotropy  or  grain  directionality  is  introduced  in  the  calendering 
process.  The  calendered  specimens  are  placed  in  a mold,  and  then  compressed.  We 
assume  that  the  identification  label  is  consistent  with  respect  to  the  direction  of  rolling, 
i.e.  if  the  label  is  parallel  to  the  direction  of  rolling,  then  it  is  consistently  parallel.  The 
dumbbell  specimens  were  cut  both  parallel  and  perpendicular  to  the  identification  label 
in  order  to  test  whether  the  tensile  response  was  sensitive  to  direction.  Figure  3.2  shows 
the  sheets. 

The  dimensions  of  the  dumbbell  specimens  are  shown  in  Figure  3.3.  These 
were  cut  using  the  Standard  Die  C (ASTM  D 412-92),  and  provided  a reduced  section 
width  of  6.35  mm  (1/4  in).  The  thickness  of  each  specimen  was  measured  in 
accordance  with  ASTM  D412-92,  three  measurements  being  made,  one  at  the  center 
and  one  at  each  end  of  the  reduced  section.  The  average  of  these  figures  was  used  as 
the  thickness.  The  specimens  were  generally  about  2.0  mm  thick.  The  difference 
between  the  maximum  and  minimum  thicknesses  was  less  than  0.08  mm  (0.003  in)  for 
each  specimen.  Consequently,  no  specimens  were  discarded,  in  accordance  with  ASTM 
D4 12-92. 


3.3. 1.1  Rubber  Uniaxial  Tension  Testing  Procedure 

Testing  was  performed  on  an  Instron  1125  test  machine  using  an  Instron  222.4  N 
(50  lb)  load  cell.  MTS  TestWorks  Software  Version  3.05-9  was  used  to  control  the 
velocity  of  the  load  head.  This  software  automatically  records  load  and  displacement, 
however,  the  recorded  displacement  is  the  grip  separation.  This  displacement 
measurement  is  undesirable  because  it  includes  end  effects  in  the  dumbbell  specimen. 

To  improve  the  strain  measurement,  gage  marks  were  placed  in  the  reduced  section, 
equidistant  from  the  center  and  perpendicular  to  the  longitudinal  axis,  and  displacement 
readings  were  manually  recorded  using  a caliper  with  a digital  measurement  display 
accurate  to  0.01  mm.  After  a displacement  and  corresponding  force  data  pair  were 
recorded,  the  caliper  opening  was  increased  close  to  the  next  desired  displacement 
reading.  In  practice,  displacement  readings  differed  from  sample  to  sample,  since  the 
important  thing  was  to  get  reasonably  spaced  measurements  with  accurate  force 
readings  at  those  displacements,  and  not  to  obtain  the  readings  at  absolutely 
predetermined  displacements.  Setting  the  caliper  at  an  exactly  precise  displacement 
was  time  consuming,  especially  for  the  tests  done  at  faster  strain  rates.  Consequently, 
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the  recorded  displacements,  hence  strains,  differed  from  one  specimen  to  another  over 
the  entire  strain  range. 

No  mechanical  conditioning  of  the  specimens  was  done.  Tests  were  performed 
at  four  strain  rates,  as  shown  in  the  Table  3.1.  In  the  actual  bearing  tests,  the 
compressive  nominal  strain  rates  varied  from  0.0003 16m/m/s  to  0.000861 m/m/s. 

Table  3.1  - Strain  rates  for  rubber  tensile  tests 


Head  Velocity  (mm/min) 

Strain  Rate  (per  second) 

2.282 

0.000499 

10.000 

0.002187 

25.400 

0.005556 

200.000 

0.043745 

3.3. 1.2  Rubber  Uniaxial  Tension  Testing  - Observations  and  Results 

Figure  3.4  shows  the  stress-strain  plots  of  all  tensile  tests.  There  is  no  clear 
difference  in  response  between  specimens  cut  parallel  to  the  identification  label  (“Hor” 
specimens)  and  those  cut  perpendicular  to  the  identification  label  (“Vert”  specimens), 
nor  by  strain  rate.  The  specimen  Batch l.Hor#2  was  chosen  because  the  stress-strain 
curve  represente  about  an  average  over  all  the  curves.  Averaging  of  all  teste  was  not 
done,  since  this  would  have  required  curve  fitting  each  test,  since  the  stress  values  of 
each  test  were  at  different  strains,  as  discussed  in  the  preceding  section. 

3.3.2  Rubber  Uniaxial  Compression  Testing 

The  B-Bearings  were  tested  to  a bearing  engineering  compressive  stress  of  just 
under  100  MPa,  so  the  uniaxial  compressive  teste  were  carried  to  a stress  level 
exceeding  100  MPa  (actually,  to  approximately  125  MPa).  At  this  stress  level,  the 
engineering  compressive  strain  exceed  0.85  m/m.  This  ensured  that  the  strain  energy 
function  would  adequately  represent  the  severe  compressive  stresses  and  strains 
expected  in  the  finite  element  analysis. 

The  bearing  manufacturer  provided  samples  normally  used  for  the  ASTM 
Compression  Set  Test  (ASTM  D 395-89).  These  samples  were  labeled  CompSet  #.  A 
second  set  were  samples  designed  to  test  the  hardness  of  the  rubber  by  use  of  the  Shore 
A durometer  (ASTM  D 2240-91),  and  were  labeled  Duro  #.  Figure  3.5  shows  the 
dimensions  and  shape  factors  of  each  type  of  specimen. 

3.3.2. 1 Rubber  Uniaxial  Compression  Testing  Procedure 

Compression  testing  was  performed  on  the  Instron  1125  using  a 98.1  KN 
(10,000  kg-f ) load  cell.  The  vertical  displacement  was  recorded  electronically,  using 
MTS  TestWorks  Software.  In  order  to  get  accurate  displacement  readings,  the  machine 
compliance  had  to  be  removed  from  the  displacement  data.  The  machine  compliance  is 
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the  measured  vertical  displacement  that  is  due  solely  to  the  flexibility  of  the  testing 
machine.  Compliance  is  significant  in  this  test  because  the  rubber  samples  resisted  a 
high  load.  Compliance  was  found  by  compressing  the  platens  together,  without  the 
specimen  in  place,  using  the  same  configuration  as  if  a specimen  were  being  tested. 

The  resulting  data  are  plotted  as  displacement  versus  load,  and  the  displacement 
(compliance)  is  found  as  a function  of  load  using  a polynomial  curvefit.  The 
compliance  is  removed  at  each  load  value  in  the  actual  test  data.  At  109  MPa,  the 
compliance  accounts  for  7.2%  of  the  measured  displacement  for  the  Compression  Set 
specimen  and  14.9%  for  the  Durometer  specimen. 

The  uniaxial  compression  test  must  simulate  homogeneous  compression.  In 
homogeneous  deformation,  straight  lines  before  deformation  remain  straight  after 
deformation  (Eringen,  1989).  In  the  compression  test  performed,  this  requires  the 
loaded  ends  to  slide  freely  along  the  loaded  surface,  so  there  is  no  bowing  of  the 
unloaded  surface.  Compression  with  the  top  and  bottom  surfaces  bonded  to  the  platens 
is  the  other  extreme  - there  is  no  sliding  along  the  loaded  surface.  The  stress-strain 
response  of  a bonded  specimen  is  a function  of  the  shape  factor.  The  larger  the  shape 
factor,  the  stiffer  is  the  response.  Figure  3.6  shows  the  difference  between 
homogeneous  and  bonded  compression. 

Attaining  perfect  homogeneous  compression  requires  a frictionless  contact 
between  the  platen  and  sample  surface.  In  order  to  simulate  frictionless  conditions  as 
closely  as  possible,  Apiezon  grease  was  applied  to  the  loaded  surfaces  prior  to  the  test. 

The  load  head  platens  had  concentric  grooves  in  them,  making  it  difficult  to 
simulate  a frictionless  surface.  Smooth  high  strength  steel  plates  were  used  between 
the  platen  and  sample  to  reduce  friction.  The  upper  plate  could  not  be  attached  to  the 
upper  load  head.  This  made  balancing  of  the  smooth  platen  difficult,  especially  since 
the  specimen  was  greased,  and  necessitated  a slight  preload  to  keep  the  upper  platen 
from  sliding  to  one  side  or  the  other. 

3.3.2.2  Rubber  Uniaxial  Compression  Testing  - Observations  and  Results 

Bowing  of  the  unloaded  surface  was  very  difficult  to  determine  visually  because 
grease  squeezing  out  from  the  edges  of  the  loaded  surfaces  obscured  the  view.  Figure 
3.7  shows  the  stress-strain  plots  of  CompSet  3 and  CompSet  4 prior  to  removal  of  the 
machine  compliance.  These  two  virtually  identical  results  show  repeatability  in  the  test, 
indicating  more  reliability  in  the  results. 

Figure  3.8  shows  the  stress-strain  curves  of  CompSet  4 and  Duro  1 before  and 
after  removal  of  the  machine  compliance.  The  "CompSet"  specimen  refers  to  a 
specimen  that  was  originally  made  for  the  ASTM  compression  set  test  (ASTM  D 395- 
89),  and  the  "Duro"  specimen  refers  to  a specimen  originally  made  for  the  ASTM 
durometer  (hardness)  test  (ASTM  D 2240-91).  The  curves  are  virtually  identical  after 
removal  of  machine  compliance,  and  indicates  a good  approximation  to  homogeneous 
compression,  since  the  shape  factors  are  different.  The  Duro  1 plot  is  shifted  more  by 
removal  of  the  compliance  because  it  is  a thinner  specimen  than  CompSet  4,  so  each 
compliance  value  is  a higher  vertical  strain.  CompSet  4 was  chosen  for  input  into  the 
finite  element  program. 
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3.3.3  Volumetric  Testing  of  Rubber 


The  volumetric  testing  was  performed  on  the  Instron  1125  using  a 98.1  KN 
(10,000  kg-f ) load  cell,  the  same  load  cell  used  for  compression  testing.  Figure  3.9 
shows  a diagram  of  the  device  that  was  manufactured  for  volumetric  testing.  The 
cylindrical  hole  was  made  at  the  specified  diameter  in  order  to  ensure  a snug  fit  with  the 
available  samples.  The  Compression  Set  samples  CompSet  5,  CompSet  6,  and 
CompSet  7 were  tested  in  this  apparatus. 

The  tight  tolerances  were  necessary  because  of  the  extreme  importance  of 
providing  no  room  for  the  rubber  to  extrude  during  the  test.  The  piston  slid  freely 
through  the  hole  when  the  bottom  plate  was  removed,  but  with  the  bottom  plate 
attached,  entrapped  air  could  not  move  around  the  piston  quickly,  thereby  retarding  its 
movement.  This  was  not  a problem,  since  the  tests  were  run  at  low  head  speeds. 

As  with  the  compression  testing,  compliance  tests  were  used  to  remove  the 
effects  of  the  inherent  flexibility  of  the  machine.  This  flexibility  was  removed  in  the 
same  fashion  as  described  above  for  the  compression  testing  (c.f.  Section  3.3.2. 1). 

Load  was  applied  to  the  piston,  and  the  load  and  corresponding  vertical 
displacement  of  the  piston  were  measured  electronically.  Figure  3.10  shows  the  load 
vs  extension  for  each  of  the  three  samples  tested  prior  to  removal  of  the  machine 
compliance.  The  slopes  are  virtually  identical.  Each  specimen  has  a very  low  initial 
stiffness,  which  appears  as  a horizontal  line.  This  initial  low  stiffness  is  probably  due  to 
the  rubber  squeezing  into  small  voids  because  the  specimens  weren’t  perfectly 
cylindrical. 

After  removal  of  the  machine  compliance,  the  raw  data  were  converted  into 
Pressure  Stress  and  Volume  Ratio.  Figure  3.11  shows  the  plot  for  CompSet  7.  Note 
that  removal  of  the  machine  compliance  does  not  account  for  the  initial  low  stiffness  in 
the  response. 

3.3.4  Rubber  Model  used  in  Finite  Element  Analysis  (based  on  Experimental 
Data) 


A particular  set  of  stress-strain  points  from  the  tensile  test  Batch l.Hor  #2  and 
the  compression  test  CompSet  4 were  chosen  for  input  into  the  finite  element  program. 
This  set  of  uniaxial  engineering  stress-strain  data  points  is  referred  to  as  Rubber  Model 
1 (RM1),  and  is  shown  in  Figure  3.12.  As  discussed  in  Section  3.2,  this  data  is  used  in 
an  ABAQUS  analysis  to  determine  the  coefficients  for  the  deviatoric  portion  of  the 
strain  energy  functions,  specifically  the  Cv  for  the  Polynomial  Strain  Energy  Function, 

and  the  ju,  and  a,  for  the  Ogden  Strain  Energy  Function. 

To  determine  the  coefficients  £>,  used  in  determining  the  form  of  the  volumetric 
portion  of  both  strain  energy  functions,  data  from  test  CompSet  7 is  used,  and  is 
referred  to  subsequently  as  CS7.  Figure  3.13  shows  the  actual  volumetric  data  used  in 
the  finite  element  model.  The  low  initial  stiffness  has  been  removed  from  the  data.  As 
discussed  in  Section  3.3.3,  the  low  initial  stiffness  is  probably  due  to  rubber  squeezing 
into  small  voids  since  the  specimens  weren’t  perfectly  cylindrical. 
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3.3.5  Measurement  Accuracy  and  Variability  and  Effect  on  Bearing  Response 

The  accuracy  of  the  Instron  1125  test  machine  is  established  through  a 
calibration  using  known  dead  loads,  and  was  calibrated  according  to  ASTM  E4,  which 

requires  that  the  machine  reading  be  within  ± 1%  of  the  true  reading.  The  machine  has 
met  these  requirements.  This  is  equivalent  to  a relative  expanded  uncertainty  (k=2)  of 
1%,  and  is  small  compared  to  specimen  variability  for  the  rubber  uniaxial  testing. 

Rubber  uniaxial  tensile  testing  was  performed  on  fourteen  (14)  dumbbell 
specimens.  The  variation  about  the  mean  of  these  tests  was  approximately  ± 10%,  so 
the  variability  can  be  attributed  to  factors  outside  the  resolution  of  the  testing  machine. 
Sources  of  variation  include  (1)  differences  in  the  properties  of  the  rubber  and  (2) 
measurement  error.  Displacements  were  recorded  manually  using  a hand-held  caliper, 
and  force  readings  were  manually  taken  at  the  time  of  each  displacement  reading. 
Measurement  error  in  the  cross-sectional  area  can  also  affect  the  results.  Based  on  this 
limited  testing  sample  and  manual  testing  methodology,  it  is  difficult  to  say  how  much 
of  the  variability  is  due  strictly  to  differences  in  the  properties  of  the  rubber,  but  it  is 
likely  that  the  largest  portion  of  the  variability  is  due  to  the  manual  measurement 
process. 

Rubber  uniaxial  compression  testing  was  performed  on  cylindrical  bullets  of 
different  shape  factors,  and  the  data  was  recorded  electronically,  as  discussed  in  Section 
3.3.2  above.  For  the  two  tests  reported  and  shown  graphically  in  Figure  3.8,  the 
variation  of  the  test  values  about  the  mean  of  the  two  tests  was  approximately  ± 2.8%. 

Volumetric  testing  was  performed  on  cylindrical  bullets,  and  the  data  was 
recorded  electronically,  as  discussed  in  Section  3.3.3  above.  For  the  three  tests  reported 
and  shown  graphically  in  Figure  3.10,  the  slopes  of  each  graph,  which  correspond  to  the 

bulk  modulus,  differ  by  only  ± 0.86%  from  the  mean  value.  This  difference  cannot  be 
attributed  to  any  specific  errors  in  testing  or  differences  in  the  bulk  modulus  from 
sample  to  sample  because  these  values  lie  inside  the  resolution  of  the  testing  machine. 

The  finite  element  analysis  of  the  bearings,  discussed  in  detail  in  Section  6, 
shows  that  the  hydrostatic  stresses  in  the  rubber  are  two  orders  of  magnitude  larger  than 
the  deviatoric  stresses.  The  uniaxial  tests  serve  to  define  the  deviatoric  strain  energy,  so 
any  variation  in  the  uniaxial  testing  results  would  have  a negligible  affect  on  the  bearing 
compressive  stress-strain  results.  The  hydrostatic  stresses  play  the  major  role  here,  but 
the  variation  in  the  three  tests  is  so  small  that  we  cannot  attribute  differences  in  the 
bearing  responses  to  differences  in  the  bulk  modulus  of  the  rubber  from  sample  to 
sample. 

3.4  Extension  of  Rubber  Experimental  Stress-Strain  Data  to  Equibiaxial  and 
Planar  Deformation  Modes  using  the  Valanis-Landel  Function 

Rubber  testing  was  performed  only  in  the  uniaxial  tensile  and  compressive 
modes.  A mode  of  deformation  independent  of  either  uniaxial  tension  or  uniaxial 
compression  is  the  planar  deformation  mode.  No  planar  deformation  tests  were  done, 
since  a test  apparatus  for  this  mode  of  deformation  was  not  available.  Planar 
deformation  is  referred  to  as  a pure  shear.  A simple  shear  is  equivalent  to  a pure  shear 
plus  a rotation.  Since  shearing  stresses  are  present  at  the  rubber/steel  interface,  accurate 
determination  of  these  stresses  is  important.  Since  it  was  not  clear  apriori  that  uniaxial 
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stress-strain  curves  would  result  in  a strain  energy  function  that  would  properly 
characterize  this  mode  of  deformation,  we  used  the  Valanis-Landel  function  to  generate 
stress  data  for  planar  deformation  as  if  we  had  performed  the  tests  ourselves.  Also, 
equibiaxial  tension,  which  is  equivalent  to  uniaxial  compression  for  an  incompressible 
material,  was  also  calculated  from  the  Valanis-Landel  function. 

The  inclusion  of  planar  tension  data  might  result  in  a strain  energy  function 
which  better  describes  the  behavior  of  the  rubber  in  all  deformation  modes.  This 
section  describes  in  detail  how  the  Valanis-Landel  function  is  found  and  used  to 
generate  these  other  data  sets.  In  Section  3.5,  the  reasonableness  of  the  planar  and 
equibiaxial  stress-strain  data  is  evaluated.  In  Section  4,  finite  element  runs  are 
presented  to  see  how  well  the  ABAQUS  code  predicts  the  experimental  and  Valanis- 
Landel  generated  data. 

For  an  elastic,  isotropic,  and  incompressible  material,  stress  differences  on  a unit 
cube  of  material  can  be  written  as 


. dU  . dU 
h hdK 

(3.4a) 

. dU  dU 

(3.4b) 

. dU  . dU 
— A, — —A, — 

^ 1 ^ dX^  4 dXl 

(3.4c) 

where  directions  1,  2,  and  3 correspond  to  a right  handed  orthogonal  coordinate  system, 
r,  is  the  true  stress  on  the  face  perpendicular  to  the  i-direction,  and  Af  is  the  stretch 
ratio,  defined  as  the  ratio  of  the  length  of  a line  in  the  deformed  configuration  to  the 
length  of  that  line  in  the  undeformed  configuration  which  initially  points  in  the  i- 
direction. 

As  noted  in  Section  3.2,  Valanis  and  Landel  (1967)  postulate  that  the  strain 
energy  function  for  rubber  can  be  written  as  a separable  function  of  the  stretch  ratios,  as 
follows: 


f/  = M(A1)  + M(A2)  + M(A3)  (3.5) 

where  w(-)  is  termed  the  Valanis-Landel  function.  Taking  the  derivative  with  respect  to 
the  stretch,  we  get 


dU  d\u{\)  1 
d\  ~ dX, 


or  u'{\) 


and  Equation  (3.4)  can  be  rewritten  as 


tj  -t2  = A1m'(A1)-  X2u\?i2) 
ti~h  = X^u'iX^ ) 


(3.6) 


(3.7a) 

(3.7b) 
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h-h  = A“'(A )-  A«'(A) 


(3.7c) 


For  a planar  tension  test  with  loading  in  the  1 direction  and  no  deformation  in 
the  2 direction,  A = 1 and  t3  = 0 . Applying  the  incompressibility  constraint 

AjAA  = 1 , we  get  A3  = — . Thus,  Equation  (3.7c)  becomes 


A \A/  (3.8) 

To  convert  Equation  (3.8)  to  engineering  stress,  note  that 

force  original  area  original  area 

= * = oi  * , where  cr  = engineering  stress. 

original  area  current  area  current  area 

In  terms  of  stretch  ratios,  tt  = cr  * — — = cr  A,  , and  the  engineering  stress  in  the  planar 

Xjkk 

mode,  in  terms  of  the  derivative  of  die  Valanis-Landel  function  , is 


1 / 1 \ 

Oi=u(Aj)-— yw  — 

A \ A / 


(3.9) 


For  an  equibiaxial  tension  test  with  loading  in  the  1 and  2 directions,  = 0 and 

tx  = t2.  In  this  case.  Equation  (3.7a)  is  identically  satisfied,  and  Equations  (3.7b)  and 
(3.7c)  give  the  result 


t2  = X2u\Xz)-X3u\X3).  (3.10) 

Applying  the  incompressibility  constraint,  and  noting  that  A = Ap  we  get 

(3.11) 

A \A  / 

Therefore,  engineering  stress  in  the  equibiaxial  mode,  in  terms  of  the  derivative  of  the 
Valanis-Landel  function,  is 


o2  -«'( A)- 


/ 

\ 


(3.12) 


Once  the  derivative  of  the  Valanis-Landel  function  is  determined,  then  the 
stress-strain  curves  for  equibiaxial  and  planar  tension  can  be  computed.  This  function 
can  be  found  from  the  existing  uniaxial  test  data.  For  a uniaxial  test  with  loading  in  the 


1 -direction,  A = A • Applying  the  incompressibility  constraint  A = A = “7=  > and 

VA 


noting  that  t2  = 0 , Equation  (3.7a)  becomes 


(3.13) 
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This  form,  however,  does  not  lend  itself  to  an  easy  identification  of  the 
derivative  of  the  Valanis-Landel  function.  A solution  to  this  problem  was  found  by 
Kearsley  and  Zappas  (1980),  as  follows: 

Rewriting  Equation  (3.13)  with  the  argument  Aj^2  yields 

h(X? ) = V2«'(^2 ) - (Ar«2)'V2«'((  A,-1'2  y* ) = X?)-  A Tu  ’(XT ) (3.14) 

Adding  Equations  (3.13)  and  (3.14)  together,  one  gets 

'iM  + '.W'1)*  A,«'(A1)-A,%'W4)  (3-15) 

Note  that  the  two  middle  terms  cancel  out.  Continuing  this  process,  we  get  two 
additional  equations, 

t,(Af ) = %\'(XT)  - X*u'(x;*),  and  (3.16) 

r,(A^)=  (3.17) 

and  adding  these  two  equations  to  Equation  (3.15),  we  get 

»i(A,)  + »1(A^)+^(A?‘)  + f1(^)-  A,«'(A1)-AfV(A?“)  (3.18) 

Grouping  the  first  and  third  terms,  and  the  second  and  fourth  terms,  we  can 
rewrite  Equation  (3.18)  as  Y ) + f1(A^(V2)(l#4)  ] = Aj*16).  We  can 

continue  this  process  and  obtain  the  series 

g[f](rO+^_(*x,'4,‘)]= v'U)-  (3.i9) 

Taking  the  limit  of  this  expression  as  n -*  oo  , Equation  (3.19)  becomes 

tog[«i(Ar* ) + 1.( X?™* )]  = \u'(\)  - «'(1)  (3.20) 

Since  A=  1 corresponds  to  no  deformation,  u'( l)  = 0 . 

The  right  hand  side  of  Equation  (3.20)  (Aw'(A))  is  an  infinite  sum  of  the  true 
stresses  at  the  stretches  specified  in  the  bracketed  expression,  but  this  sum  converges 
rapidly.  Note  that  the  first  and  second  terms  are  on  opposite  sides  of  the  undeformed 

state  (A  =1),  so  that  both  extension  and  compression  data  are  required,  and  the  series 
terms  oscillate  around  the  undeformed  state  and  rapidly  approach  it 

In  calculating  Equation  (3.20),  the  same  uniaxial  data  set  was  used  as  for  RM1; 
tension  data  from  Batch l.Hor#2  and  compression  data  from  CompSet  4.  Every  data 
point  in  this  set  was  used.  The  experimental  uniaxial  data  is  converted  to  true  stress- 
strain  data  by  multiplying  each  stress  by  its  corresponding  stretch.  Figure  3.14  shows 
the  true  stress-strain  plot.  To  aid  in  the  calculations,  it  was  necessary  to  obtain  a 
functional  representation  of  the  uniaxial  stress-strain  curve.  Therefore,  the  curve  was 
fitted  with  polynomial  functions  found  from  a least-square  curvefitting  procedure.  Due 
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to  the  sharp  curvature,  it  was  not  possible  to  adequately  fit  the  curve  with  one 
polynomial  function  over  its  entire  domain.  Instead,  the  domain  was  broken  up  into 
segments,  each  of  which  was  curvefitted  separately.  Once  completed,  a computer 
program  was  written  to  find  the  values  u'(k)  at  desired  stretch  values. 


Figure  3.15  shows  the  polynomial  curvefit  over  a particular  range  of  stretch 
values.  A fifth  order  polynomial  was  used  over  this  range.  Over  some  other  ranges,  a 

ninth  order  polynomial  was  used.  Figure  3.16  is  a plot  of  w'(A) . A listing  of  the 
program  used  to  find  these  values  is  shown  in  Appendix  B.  I 


Table  3.2  is  the  output  from  a particular  run  of  the  program  used  to  calculate 
w'(A)  and  shows  that  convergence  within  0.0001  occurs  in  less  than  11  iterations. 


Table  3.2  - Calculated  values  of  the  derivative  of  the  Valanis-Landel  Function  w'(A) 
Tolerance  = 0.000100 


Max 

#Iters 

Stretch 

Sum[i] 

Sum[i-1] 

Difference 

100 

10 

0.15 

-67.7246 

-67.7246 

0.0000 

100 

9 

0.35 

-0.3872 

-0.3871 

-0.0001 

100 

8 

0.55 

0.2652 

0.2653 

-0.0001 

100 

8 

0.75 

0.1148 

0.1148 

-0.0000 

100 

6 

0.95 

0.0048 

0.0049 

-0.0001 

100 

7 

1.15 

0.7627 

0.7627 

0.0001 

100 

8 

1.35 

1.0948 

1.0948 

0.0000 

100 

8 

1.55 

1.3827 

1.3826 

0.0000 

100 

8 

1.75 

1.6109 

1.6109 

0.0000 

100 

8 

1.95 

1.9006 

1.9005 

0.0000 

100 

8 

2.15 

2.2284 

2.2284 

0.0000 

100 

8 

2.35 

2.5715 

2.5714 

0.0000 

100 

8 

2.55 

2.9837 

2.9837 

0.0000 

100 

8 

2.75 

3.4145 

3.4145 

0.0000 

100 

8 

2.95 

3.8874 

3.8874 

0.0000 

100 

8 

3.15 

4.3750 

4.3749 

0.0000 

100 

8 

3.35 

4.8754 

4.8754 

0.0000 

100 

8 

3.55 

5.3961 

5.3961 

0.0000 

100 

8 

3.75 

5.9439 

5.9439 

0.0000 

100 

8 

3.95 

6.5195 

6.5195 

0.0000 

100 

8 

4.15 

7.1172 

7.1171 

0.0000 

100 

8 

4.35 

7.7281 

7.7281 

0.0000 

100 

8 

4.55 

8.3451 

8.3450 

0.0000 

100 

8 

4.75 

8.9669 

8.9668 

0.0000 

100 

8 

4.95 

9.5974 

9.5974 

0.0000 

100 

8 

5.15 

10.2452 

10.2452 

0.0000 

100 

8 

5.35 

10.9185 

10.9185 

0.0000 

100 

8 

5.55 

11.6196 

11.6196 

0.0000 

100 

8 

5.75 

12.3393 

12.3393 

0.0000 

100 

8 

5.95 

13.0581 

13.0581 

0.0000 

100 

8 

6.15 

13.7550 

13.7549 

0.0000 

100 

8 

6.35 

14.4305 

14.4305 

0.0000 

40 


3.5  Assessment  of  the  Valanis-Landel  Function 

Equation  (3.13)  gives  the  true  stress  for  a uniaxial  test.  Converting  true  stress 
to  engineering  stress,  we  get 


<7  = u\X)  - A"3/V(r,/2).  (3.21) 

Figures  3.17  and  3.18  show  the  comparison  between  the  experimental  data  and 
the  engineering  stresses  generated  by  Equation  (3.21)  for  uniaxial  tension  and 
compression,  respectively.  The  results  are  in  agreement  with  the  experimental  data, 
giving  us  confidence  in  the  evaluation  of  the  derivative  of  the  Valanis-Landel 
function. 


The  stress-strain  data  for  equibiaxial  tension  were  generated  from  Equation 
(3.12).  Since  uniaxial  compression  and  equibiaxial  tension  are  equivalent  for  an 
incompressible  material,  we  can  transform  the  experimental  uniaxial  compression  data 
to  equibiaxial  tension  and  compare  the  results  to  the  equibiaxial  stresses  predicted  by 
the  Valanis-Landel  function.  Figure  3.19  shows  the  equivalence  between  uniaxial 
compression  and  equibiaxial  tension  (ABAQUS  1995a).  The  stresses  in  Figure  3.19 
are  true  (Cauchy)  stresses,  and  the  superposition  of  a hydrostatic  tensile  stress  -t  on  the 
uniaxial  compressive  stress  t results  in  an  equibiaxial  tensile  stress  -t,  without  any 
change  in  deformation. 


With  reference  to  Figure  3.19,  the  incompressibility  assumption  gives  us,  in 
terms  of  the  stretches,  the  equation  AXAVAZ  = 1 . Using  this  relationship,  the  stretch  in 

the  coordinate  directions  x and  y is  A = jy- , and  the  equibiaxial  engineering  strain  is 

\*z 


e = A- 1 = 


(3.22) 


The  true  uniaxial  compressive  stress  t,  in  terms  of  the  engineering  uniaxial 
compressive  stress  <7Z,  is 


t = (3.23) 

and  the  engineering  equibiaxial  stress,  in  terms  of  the  true  stress,  is 

ax  = ay\.  (3.24) 

Noting  that  ax  = <7y,  and  denoting  this  stress  as  cr equibiaxial,  by  combining 
Equations  (3.22),  (3.23),  and  (3.24), 

G equibiaxial  = & z^z  ' (3.25) 

Equation  (3.25)  represents  the  equibiaxial  stress-strain  response  that  is  based 
on  the  assumptions  that  the  rubber  is  isotropic  and  incompressible.  Figure  3.20  shows 
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excellent  agreement  between  the  equibiaxial  stresses  calculated  from  Equation  (3.25) 
and  the  stresses  found  from  the  Valanis-Landel  function. 

The  stress-strain  data  for  planar  tension  were  generated  from  Equation  (3.9). 
Figure  3.21  is  a comparison  between  the  planar  tensile  stresses  and  the  experimental 
uniaxial  tensile  stresses,  and  it  is  seen  that  the  planar  stresses  are  at  least  as  large  as  the 
uniaxial  stresses.  The  relative  values  are  consistent  with  experimental  data  reported 
by  Treloar  (1943),  as  shown  in  Figure  3.22.  Furthermore,  the  force  required  to  deform 
a specimen  in  planar  tension  should  be  at  least  as  large  as  the  force  required  in 
uniaxial  tension  because  in  the  planar  deformation  mode,  the  specimen  is  fixed  along 
an  axis  perpendicular  to  the  direction  of  loading,  whereas  in  the  uniaxial  mode,  the 
specimen  is  free  to  move  in  this  same  direction. 

3.6  Rubber  Models  used  in  Finite  Element  Analysis  (based  on  Experimental  and 
Valanis-Landel  Generated  Data) 

Figure  3.23  shows  the  stress-strain  curve  for  planar  tension,  as  found  from  the 
Valanis-Landel  function,  together  with  the  experimental  uniaxial  tensile  and 
compressive  stress-strain  curve  (RM1).  Each  of  these  curves  represents  an 
independent  mode  of  deformation.  Collectively,  this  stress-strain  data  set  is  termed 
the  VL1  rubber  model. 


A third  rubber  model,  termed  VL2,  is  shown  in  Figure  3.24.  This  model 
differs  from  VL1  by  substituting  the  equibiaxial  data  generated  from  the  Valanis- 
Landel  function  for  the  experimental  uniaxial  compression  data. 


Stresses  for  a stretch  greater  than  2.64  (strain  of  1.64)  for  equibiaxial  tension 
cannot  be  found  because  the  lowest  stretch  value  in  the  uniaxial  tests  is  0.14389. 
Higher  stretches  exceed  the  limits  of  the  underlying  data  set,  i.e.  from  the  second  term 

1 ^ f— i^-l  = M'(0.143). 

U-64 2J  v ' 


of  Equation  (3.12);  u'\  -j  ] = u 


Each  rubber  model,  VL1  and  VL2,  is  used  by  ABAQUS  to  determine  the 
coefficients  for  the  deviatoric  portion  of  the  strain  energy  functions,  specifically  the  Qj 
coefficients  for  the  Polynomial  Strain  Energy  Function,  and  the  /i,  and  a,  coefficients 
for  the  Ogden  Strain  Energy  Function. 

In  Section  4. 1 , finite  element  runs  are  presented  to  show  how  well  the  finite 
element  program  predicts  the  experimental  and  Valanis-Landel  generated  stress-strain 
curves,  using  the  deviatoric  rubber  models  RM1,  VL1,  and  VL2,  and  to  show  how 
well  the  finite  element  program  predicts  the  experimental  volumetric  compression  test 
through  use  of  volumetric  rubber  model  CS7. 
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80  160  240  320  400 

Extension  (mm) 

Figure  3.1  - Effect  of  mechanical  conditioning  on  the  stress-strain  response  of  rubber 


Specimens  labelled  Batch#. Hor  Specimens  labelled  Batch#.Vert 

Figure  3.2  - Rubber  sheets  from  which  tensile  specimens  were  cut 
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Figure  3.3  - Dimensions  of  rubber  tensile  specimens 


▲ Batch  1. Vert  #3  - strain  rate  0.043745/sec 

□ Batch  1 .Hor#4  - strain  rate  0.043745/sec 

o Batch  l.Vert  #1  - strain  rate  0.005556/sec 

v Batch  1 .Vert  #2  - strain  rate  0.005556/sec 

o Batch  4. Vert  #1  - strain  rate  0.002187/sec 

a Batch  4.Hor#l  - strain  rate  0.002187/sec 

x Batch  2.Vert#l  - strain  rate  0.002187/sec 

* Batch  2.Hor#l  - strain  rate  0.002187/sec 

♦ Batch  3.Vert#l  - strain  rate  0.002187/sec 
Batch  3.Hor#l  - strain  rate  0.002187/sec 

— * — Composite  Batch.Vert#l-  strain  rate  0.002187/sec 

— b-  - Composite  Batch.Vert  #2-3  strain  rates 

- - * - - Batch  l.Vert  #4  - strain  rate  0.000499/sec 

Batch  1 .Hor  #2  - strain  rate  0.000499/sec 


0 1 2 3 4 5 6 7 


Figure  3.4  - Stress-strain  curves  for  all  rubber  tensile  tests 
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Engineering  Tensile  Stress  (psi) 


Figure  3.5  - Dimensions  of  rubber  specimens  used  for  compression  testing 
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Homogeneous  Compression 


Bonded  Compression 


Figure  3.6  - Homogeneous  vs  bonded  compression 
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▲ CompSet4  - Machine  Compliance  Not  Removed 
o CompSet3  - Machine  Compliance  Not  Removed 
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Figure  3.7  - Comparison  of  compressive  response  of  two  compression  set  specimens 


▲ CompSet4  - Machine  Compliance  Not  Removed 
• Durol- Machine  Compliance  Not  Removed 
A CompSet4  - Machine  Compliance  Removed 
o Durol- Machine  Compliance  Removed 
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Figure  3.8  - Comparison  of  compressive  response  of  compression  set  and  durometer 
specimens 
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Extension  data  has  been  adjusted  to  make 
zero  extension  correspond  to  zero  load, 
however,  the  machine  compliance  has  not 
been  removed  from  this  data. 


0 0.01  0.02  0.03  0.04  0.05  0.06  0.07  0.08 

Adj  Extension  (in) 

Figure  3.10  - Comparison  of  volumetric  compression  tests 


- Specimen  CompSet  5 
i Cs — Specimen  CompSet  6 
-e — Specimen  CompSet  7 


0.93  0.94  0.95  0.96  0.97  0.98  0.99  1 


Figure  3.1 1 - CompSet  7 (CS7)  volumetric  compression  test 
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Experimental  Pressure  (MPa) 
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Figure  3.12  - Rubber  model  1 (RM1) 
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Figure  3.13  - Volumetric  compression  plot  (based  on  CS7) 
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Figure  3. 14  - True  stress-strain  data  set  used  to  find  the  derivative  of  the  Valanis- 
Landel  Function  («'(A)) 
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True  Engineering  Uniaxial  Stress 


Y = M0  + M1*x  + ...  M8*x8  + M9*x9 

M0 

-1.1664817e+02 

Ml 

9.8062739e+02 

M2 

-4.5649689e+02 

M3 

-2.1791899e+04 

M4 

8.4299544e+04 

M5 

-9.6543082e+04 

R 

9.9997025e-01 

0.1  0.15  02  0.25  03  0.35 


Stretch 

Figure  3.15  - Polynomial  curvefit  used  to  obtain  a functional  representation  of  the 
stress-strain  curve  for  determination  of  the  derivative  of  the  Valanis- 

Landel  Function  (w'(A))  for  0. 14  < A < 0.3 
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Figure  3.16  - Derivative  of  the  Valanis-Landel  Function  («'(A)) 


▲ Experimental  Tensile  Stress:  Batch  l.Hor#2 
Uniaxial  Tensile  Stress  from  VL  function 
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Figure  3.17  - Comparison  between  the  experimental  uniaxial  tension  test  and  the 
uniaxial  tension  derived  from  the  Valanis-Landel  function 
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x Experimental  Compressive  Stress  - CompSet  4 
— Uniaxial  Compressive  Stress  from  VL  function 

-1  -0.8  -0.6  -0.4  -0.2  0 
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Engineering  Strain 

Figure  3.18-  Comparison  between  the  experimental  uniaxial  compression  test  and  the 
uniaxial  compression  derived  from  the  Valanis-Landel  function 


Figure  3.19  - Equivalence  of  uniaxial  compression  and  equibiaxial  tension  for  an 
incompressible  material 
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O Uniaxial  Compression  converted  to  Equibiaxial  Tension 
Equibiaxial  Tensile  Stress  found  from  VL  function 


0 0.5  1 1.5  2 


Figure  3.20  - Comparison  between  the  equibiaxial  tension  converted  from  uniaxial 
compression  and  the  equibiaxial  tension  derived  from  the  Valanis- 
Landel  function 


A Experimental  Uniaxial  Tensile  Stress:  Batch  1 .Hor  #2 
Planar  Tensile  Stress  found  from  VL  function 
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Figure  3.21  - Comparison  between  the  experimental  uniaxial  tension  test  and  the 
planar  tension  derived  from  the  Valanis-Landel  function 
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Engineering  Stress  (MPa)  Engineering  Stress  (MPa) 


▲ Uniaxial  Tension  Test  Data 
o Planar  Tension  Test  Data 
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Figure  3.22  - Treloar’s  experimental  data 


o Experimental  Uniaxial  Stress  - RM1  (MPa) 

■ Planar  Tensile  Stress  found  from  VL  function 
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Figure  3.23  - VL1  rubber  model 
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▲ Experimental  Uniaxial  Tension  : Batchl.Hor#2 
x Planar  Stress  found  from  VL  function 
o Equibiaxial  Stress  found  from  VL  function 
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Figure  3.24  - VL2  rubber  model 
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4.  INVESTIGATION  OF  THE  RUBBER  MODELS  USED  IN  THE  FINITE 
ELEMENT  ANALYSIS 


4.1  Calibration  of  the  Rubber  Models  in  ABAQUS 

To  validate  the  finite  element  modeling  of  rubber,  we  need  to  see  how  well  the 
finite  element  results  approximate  the  experimental  and  Valanis-Landel  generated 
stress-strain  curves.  In  this  section  we  select  the  number  of  terms  to  be  included  in 
each  strain  energy  function,  and  the  element  type  to  be  used  in  the  finite  element 
analysis.  The  selected  element  has  been  subjected  to  the  proper  boundary  conditions 
to  reproduce  each  of  the  four  modes  we  are  analyzing,  namely  uniaxial  tension, 
uniaxial  compression,  equibiaxial  tension,  and  planar  tension.  At  the  end  of  this 
section,  the  finite  element  results  are  compared  directly  to  the  experimental  data  and 
the  data  generated  by  the  Valanis-Landel  function  for  equibiaxial  and  planar  tension. 

To  determine  how  many  terms  from  the  strain  energy  potentials  to  include  in 
the  analysis,  it  is  helpful  to  look  at  the  reduced  stress.  For  deformations  in  uniaxial 
tension  and  compression,  the  true  stress  is 


f-2  A 


-i) 


dU 


Dividing  Equation  (4.1)  through  by  A2  - 


1 


\_dlf 
A dI2J 

defines  a measure  called  the 


"reduced"  stress  (McKenna  and  Zapas,  1983;  Yeoh,  1993),  shown  as 


(4.1) 


(a2-i/a) 


= 2 


'dU_  \_dlT 
K dlx  A dl2 1 


(4.2) 


The  terms  on  the  right-hand  side  provide  a sensitive  measure  of  the  deviations 
of  the  elastic  properties  from  constant  values  (McKenna  and  Zapas  1983).  By  plotting 

— as  the  abscissa,  the  compression  information  is  shown  in  the  region  greater  than  1, 
A 

and  the  tension  information  is  compressed  between  0 and  1. 

For  the  Polynomial  Strain  Energy  Potential,  N=1  is  the  Mooney-Rivlin  model. 
Written  out  explicitly,  we  get 


U = C10  (/j  - 3)  + C01  (/2  - 3)  (4.3) 

Taking  the  appropriate  derivatives  of  Equation  (4.3)  and  inserting  them  into 
Equation  (4.2),  the  reduced  stress-strain  relation  for  uniaxial  tension  and  compression 
is 


Ir  — 2C10  + ^ 2C01 


(4.4) 
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Equation  (4.4)  shows  that  the  Mooney-Rivlin  model  will  produce  reduced 

stresses  that  are  linear  when  plotted  against  . Figure  4. 1 shows  the  experimental 

/i 

uniaxial  tensile  and  compressive  stresses  expressed  as  reduced  stresses  and  plotted 

versus  y . It  is  clear  that  a linear  approximation  will  not  be  sufficient  to  properly 

model  the  rubber.  For  this  reason,  the  calibration  runs  were  performed  using  N-2  for 
the  Polynomial  Strain  Energy  Potential. 


This  type  of  analysis  is  not  so  easily  done  with  the  Ogden  Strain  Energy 
Potential.  The  Ogden  model  is  written  in  terms  of  the  stretches,  not  the  strain 
invariants.  However,  the  larger  the  number  of  terms  included  (N),  the  easier  the 
model  should  capture  the  nonlinearity  in  the  stress-strain  curve.  Due  to  the  high 
nonlinearity  present  in  the  experimental  data,  the  calibration  runs  used  N-6,  which  is 
the  largest  number  of  terms  available  in  the  ABAQUS  code  . Lower  values  of  N 
require  less  computational  effort  to  obtain  the  stresses,  so,  for  comparison,  the  Ogden 
model  with  N=3  was  tested.  In  Section  4.2,  the  calibration  runs  for  the  uniaxial 
deformation  mode  will  be  shown  in  terms  of  reduced  stresses. 


The  uniaxial,  equibiaxial,  and  planar  calibration  runs  were  conducted  using 
one  3-dimensional,  linear  displacement  “hybrid”  finite  element  with  reduced 
integration  (ABAQUS  element  type  C3D8RH),  with  unit  dimensions  in  X,  Y,  and  Z. 
Elements  in  ABAQUS  that  use  the  mixed  formulation  (c.f.  Section  3.2)  are  called 
hybrid  elements.  Proper  application  of  boundary  conditions  allows  us  to  perfectly 
simulate  these  homogeneous  deformations.  The  use  of  one  finite  element  with  an 
assumed  linear  displacement  is  sufficient,  since  homogeneous  deformations  have  a 
constant  strain  field  and  a single  element  with  an  assumed  linear  displacement  can 
represent  this  state  of  strain. 

For  each  of  the  four  deformation  modes,  compressible  and  incompressible 
analyses  were  carried  out.  Table  4.1  shows  the  calibration  runs  performed.  CompSet 
7 (CS7  - see  Figure  3.13)  was  the  volumetric  compression  data  used  for  all  the 
Compressibility  runs.  RM1  is  shown  in  Figure  3.12,  the  VL1  rubber  model  in  Figure 
3.23,  and  the  VL2  rubber  model  in  Figure  3.24. 
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Table  4.1  - Rubber  calibrations 


Type  of  Test 

Compressible  Models 
(uses  CS7) 

Incompressible  Models 
(excludes  CS7) 

Uniaxial  Tension 

RM1,  VL1,  VL2 

RM1 

Uniaxial  Compression 

RM1,  VL1,  VL2 

RM1 

Equibiaxial  Tension 

RM1,  VL1,  VL2 

RM1 

Planar  Tension 

RM1,  VL1,  VL2 

RM1 

Volumetric  Compression 

RM1,  VL1 

Not  Applicable 

Note:  All  tests  have  been  conducted  using  the  Ogden  Strain  Energy  Potential,  N=6 
and  N=3,  and  the  Polynomial  Strain  Energy  Potential,  N=2,  except  the 
Volumetric  Test,  in  which  only  the  Ogden  form,  N=6  was  run. 


For  the  volumetric  compression  calibration  runs,  an  axisymmetric  finite 
element  model  was  used.  This  model  is  appropriate  since  the  volumetric  specimen  is 
circular  and  the  loading  is  uniformly  distributed.  A slice  of  the  specimen  in  the  R-Z 
plane  is  modeled,  from  R = 0 to  the  outer  edge  and  from  the  midheight  to  the  top,  so 
that  1/4  of  a slice  through  the  axis  of  symmetry  is  modeled.  CAX4H  elements  were 
used,  which  are  hybrid  2-dimensional  linear  displacement  elements.  Fixed  boundary 
conditions  were  set  along  the  midheight  nodes  in  both  the  R and  Z directions,  and 
along  the  outer  edge  and  the  top  nodes  in  the  R direction.  The  axisymmetric  element 
automatically  applies  fixed  boundary  conditions  along  R = 0.  Displacement  boundary 
conditions  were  placed  on  the  top  nodes,  moving  each  node  down  by  the  same  amount 
(the  negative  Z direction),  and  the  RIKS  arc-length  method  was  used.  The  model 
converged  rapidly. 

As  shown  in  Table  4.1,  for  the  uniaxial,  equibiaxial,  and  planar  modes  of 
deformation,  the  calibrations  were  done  with  and  without  compressibility,  for  rubber 
model  RM1.  The  finite  element  results  are  virtually  identical  whether  or  not 
compressibility  is  included,  so  the  VL1  and  VL2  calibrations  were  run  using  only  a 
compressible  model.  The  similarity  between  compressible  and  incompressible  results 
is  because  the  homogeneous  deformations  allow  the  rubber  to  freely  expand  or 
contract,  so  as  to  nearly  retain  its  volume.  Hence  the  volumetric  deformation  is 
negligibly  small. 

The  RM1,  VL1,  and  VL2  rubber  models  are  used  to  determine  the  deviatoric 
parameters.  For  purposes  of  comparing  the  results  by  model,  it  is  important  to  keep  in 
mind  which  deformation  modes  are  included  in  each  of  the  models.  Table  4.2 
summarizes  the  test  data  included  in  each  model.  Following  this  table  are  average 
Root  Mean  Square  (RMS)  errors  for  each  calibration  run  and  plots  of  the  results, 
organized  by  deformation  mode. 
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Table  4.2  - Test  data  included  in  each  rubber  model 


Deformation  Mode 

RM1 

Rubber  Models 

VL1 

VL2 

Uniaxial  Tension 

X 

X 

X 

Uniaxial  Compression 

X 

X 

Equibiaxial  Tension 

X 

Planar  Tension 

X 

X 

The  average  RMS  error  is  calculated  from  the  formula  shown  in  Equation 
(4.5),  and  measures  the  relative  error  in  stress. 


avg 


RMS  error  = 


1 

N 

y 

(o  — (j  > 

Uexp  U FE 

2 

Zj 

1=1 

k ®exP  y 

(4.5) 


M is  the  number  of  points  at  which  the  stress  was  calculated  by  the  finite 
element  analysis,  aFE  is  the  stress  calculated  by  the  finite  element  analysis,  and  <Jexp  is 
the  experimentally  determined  stress. 

Table  4.3  shows  the  average  RMS  error  for  each  of  the  calibration  runs  for 
Uniaxial  Tension.  The  performance  of  each  model  is  good,  with  slightly  better  results 
from  rubber  models  VL1  and  VL2.  Figures  4.2  through  4.5  are  plots  of  the  finite 
element  results  versus  the  experimental  tension  data,  and  graphically  depict  the 
closeness  of  the  finite  element  approximations  to  the  tensile  data,  using  both 
experimental  data  and  data  generated  by  the  Valanis-Landel  function. 


Table  4.3  - Average  RMS  error  of  calibration  runs  for  uniaxial  tension 


Rubber  Model 

Ogden  (N=6) 

Ogden  (N=3) 

Poly.  (N=2) 

RM1 

Compressible 

0.0195 

0.0159 

0.0177 

RM1 

Incompressible 

0.0195 

0.0159 

0.0177 

VL1 

Compressible 

0.0119 

0.0160 

0.0195 

VL2 

Compressible 

0.0117 

0.0162 

0.0156 

Table  4.4  shows  the  average  RMS  error  for  each  of  the  calibration  runs  for 
Uniaxial  Compression.  Each  finite  element  result  provides  a good  approximation  to 
the  experimental  compression  data,  with  RM1  slightly  better.  Also,  note  that  the 
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Ogden  model  with  6 terms  (N=6)  results  in  the  lowest  average  RMS  error  for  each 
rubber  model.  Figures  4.6  through  4.9  are  plots  of  the  finite  element  results  versus  the 
experimental  compression  data. 


Table  4.4  - Average  RMS  error  of  calibration  runs  for  uniaxial  compression 


Rubber  Model 

Ogden  (N=6) 

Ogden  (N=3) 

Poly.  (N=2) 

RM1 

Compressible 

0.0280 

0.0446 

0.0476 

RM1 

Incompressible 

0.0280 

0.0446 

0.0476 

VL1 

Compressible 

0.0305 

0.0451 

0.0508 

VL2 

Compressible 

0.0309 

0.0452 

0.0417 

Table  4.5  shows  the  average  RMS  error  for  each  of  the  calibration  runs  for 
equibiaxial  tension,  and  Figures  4.10  through  4.13  show  plots  of  each  run.  The 
comparisons  are  between  the  data  generated  by  the  Valanis-Landel  function  and  the 
ABAQUS  results  using  the  three  rubber  models.  The  results  for  each  rubber  model 
are  comparable,  but  experimental  rubber  model  RM1  does  a slightly  better  job  than 
VL1  and  VL2  when  using  the  Ogden,  N=6  and  Polynomial,  N=2  strain  energy 
functions,  with  the  Ogden,  N=6  providing  the  best  fit  to  the  equibiaxial  data  generated 
by  the  Valanis-Landel  function.  It  is  important  to  note  that  rubber  models  RM1  and 
VL1  do  not  include  the  equibiaxial  data  set,  yet  ABAQUS  closely  approximates  the 
equibiaxial  deformation  mode. 

Table  4.5  - Average  RMS  error  of  calibration  runs  for  equibiaxial  tension 


Rubber  Model 

Ogden  (N=6) 

Ogden  (N=3) 

Poly.  (N=2) 

RM1 

Compressible 

0.0117 

0.0271 

0.0226 

RM1 

Incompressible 

0.0118 

0.0268 

0.0225 

VL1 

Compressible 

0.0135 

0.0273 

0.0389 

VL2 

Compressible 

0.0133 

0.0262 

0.0256 

Table  4.6  shows  the  average  RMS  error  for  each  of  the  calibration  runs  for 
planar  tension,  and  Figures  4. 14  through  4. 17  show  plots  of  each  run.  An  important 
observation  is  that  the  Ogden  strain  energy  function,  using  rubber  model  RM1  (which 
does  not  include  the  planar  data),  does  a much  better  job  of  fitting  the  planar  data  than 
does  the  Polynomial  strain  energy  function.  This  is  probably  due  to  the  fact  that  the 
Ogden  model  is  a special  form  of  the  Valanis-Landel  function.  While  this  alone  does 
not  guarantee  that  the  planar  data  are  correct,  the  planar  tensile  stresses  from  the 
Valanis-Landel  function  are  slightly  larger  than  the  experimental  uniaxial  stresses,  and 
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this  is  similar  to  the  relationship  between  the  experimentally  determined  planar 
tension  data  and  uniaxial  tension  data  reported  by  Treloar  (1943).  A second 
observation  is  that  when  the  planar  data  is  included  in  the  data  set  used  to  determine 
the  strain  energy  function  (VL1  and  VL2),  the  Polynomial  strain  energy  function 
performs  nearly  as  well  as  the  Ogden  function. 


Table  4.6  - Average  RMS  error  of  calibration  runs  for  planar  tension 


Rubber  Model 

Ogden  (N=6) 

Ogden  (N=3) 

Poly.  (N=2) 

RM1 

Compressible 

0.0182 

0.0147 

0.0954 

RM1 

Incompressible 

0.0181 

0.0144 

0.1012 

VL1 

Compressible 

0.0100 

0.0123 

0.0243 

VL2 

Compressible 

0.0097 

0.0121 

0.0165 

The  volumetric  test  runs  were  only  performed  using  the  compressibility  data 
(CS7)  for  rubber  models  RM1  and  VL1.  Using  an  incompressible  model  would  result 
in  no  deformation,  since  the  rubber  is  perfectly  confined.  As  shown  in  Figure  4.18, 
the  finite  element  program  accurately  reproduces  the  experimental  data,  and  as 
expected,  there  is  no  difference  between  the  use  of  either  RM1  or  VL1,  because  those 
models  determine  the  deviatoric  parameters,  whereas  data  from  CS7  is  used  to 
determine  the  volumetric  parameters  and  is  the  same  for  both  calibration  runs. 

4.2  Reduced  Stress  Plots 

As  noted  in  Section  4. 1 , the  reduced  stress  provides  a sensitive  measure  of  the 
deviations  of  the  elastic  properties  from  constant  values.  Reduced  stress  plots,  also 
referred  to  as  Mooney-Rivlin  plots,  are  frequently  used  by  researchers  who  study  the 
mechanical  properties  of  rubber.  The  formula  for  reduced  uniaxial  stress  is  repeated 
here  for  convenience. 


t 

(A2  - 1 /A) 


= 2 


A 


■ + 


A dl2  J 


(4.2) 


By  plotting  — as  the  abscissa,  the  compression  information  is  shown  in  the  region 

A 

greater  than  1,  and  the  tension  information  is  compressed  between  0 and  1. 


In  this  section,  the  calibration  results  for  uniaxial  stress  are  presented  in  terms 
of  the  reduced  stresses.  Figures  4.19  through  4.21  show  plots  of  the  reduced  stress 

versus  — for  the  RM1,  VL,  and  VL2  rubber  models,  respectively.  These  plots  mirror 

A 

the  results  shown  in  the  calibration  runs  in  Section  4. 1 . 
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The  experimental  data  exhibit  a rapid  increase  in  modulus  as  the  unstrained 
state  is  approached.  This  has  been  observed  by  McKenna  and  Zapas  (1983)  and  other 
authors.  Finite  element  results  do  not  capture  this  change  in  reduced  stress.  McKenna 
and  Zapas  reported  that  the  compression  modulus  is  greater  than  the  tension  modulus 
near  the  undeformed  state.  This  observation  agrees  with  our  experimental  data.  The 
reasons  for  this  phenomenon  are  not  clear.  McKenna  and  Zapas  (1983)  used  bonded 
cylinders  and  corrected  the  results  to  homogeneous  compression  using  factors 
determined  in  a finite  difference  analysis.  They  conjectured  that  this  phenomenon 
might  be  due  to  their  experimental  approach.  However,  in  our  experiments,  the  tests 
were  carried  out  on  unbonded  cylinders  that  had  been  greased  to  reduce  friction  along 
the  loaded  surfaces. 

It  should  be  noted  that  most  work  in  rubber  using  the  VL  function  and 
Mooney-Rivlin  plots  has  been  done  with  unfilled  rubber.  However,  the  rubber  used  in 
this  testing  program  is  filled  with  carbon  black.  This  filler  was  added  primarily  to 
control  the  shear  stiffness  and  damping  characteristics.  The  Mooney-Rivlin  plots  of 
unfilled  rubbers  exhibit  different  behavior  than  is  observed  here.  With  unfilled 
rubbers,  the  response  starts  at  a low  value,  goes  toward  a maximum  in  the  vicinity  of 

= 1.25  and  then  decreases  mildly  before,  possibly,  increasing  again.  Here  the 
response  at  values  of  Y\  > 1-25  is  strongly  concave  upwards. 
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Reduced  Tensile  Stress  (MPa),  Experimental 
Reduced  Compressive  Stress  (MPa),  Experimenta 


0 1 2 3 4 5 6 7 


Figure  4. 1 - Plot  of  reduced  stresses  for  the  experimental  uniaxial  tensile  and 
compressive  stress-strain  data 
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Figure  4.2  - Calibration  to  uniaxial  tension  with  compressibility  using  RM  1 
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Figure  4.3  - Calibration  to  uniaxial  tension  without  compressibility  using  RM1 
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Engineering  Tensile  Stress  (MPa)  4^  Engineering  Tensile  Stress  (MPa) 
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4 - Calibration  to  uniaxial  tension  with  compressibility  using  VL1 
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Figure  4.5  - Calibration  to  uniaxial  tension  with  compressibility  using  VL2 
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Figure  4.6  - Calibration  to  uniaxial  compression  with  compressibility  using  RM 1 
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Figure  4.7  - Calibration  to  uniaxial  compression  without  compressibility  using  RM1 
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Figure  4.9  - Calibration  to  uniaxial  compression  with  compressibility  using  VL2 
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Equibiaxial  Stress  found  from  VL  function 

O Equibiaxial  Tension,  Og,  N=6,  RM1 , Compressible 
A Equibiaxial  Tension,  Og,  N=3,  RM 1 , Compressible 
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Figure  4.10  - Calibration  to  equibiaxial  tension  with  compressibility  using  RM1 
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o Equibiaxial  Tension,  Og,  N=6,  RM  1 , Incompressible 
A Equibiaxial  Tension,  Og,  N=3,  RM  1 , Incompressible 
x Equibiaxial  Tension,  Po,  N=2,  RM  1 , Incompressible 
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Figure  4.1 1 - Calibration  to  equibiaxial  tension  without  compressibility  using  RM1 
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Equibiaxial  Stress  found  from  VL  function 

o Equibiaxial  Tension,  Og,  N=6,  VL1,  Compressible 
A Equibiaxial  Tension,  Og,  N=3,  VL1,  Compressible 
x Equibiaxial  Tension,  Po,  N=2,  VL1 , Compressible 
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Figure  4.12  - Calibration  to  equibiaxial  tension  with  compressibility  using  VL1 
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Figure  4.13  - Calibration  to  equibiaxial  tension  with  compressibility  using  VL2 
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Planar  Stress  found  from  VL  function 

o Planar  Tension,  Og,  N=6,  RM 1 , Compressible 
A Planar  Tension,  Og,  N=3,  RM  1 , Compressible 
x Planar  Tension,  Po,  N=2,  RM  1 , Compressible 
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Figure  4.14  - Calibration  to  planar  tension  with  compressibility  using  RM1 
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Planar  Stress  found  from  VL  function 

O Planar  Tension,  Og,  N=6,  RM  1 , Incompressible 
A Planar  Tension,  Og,  N=3,  RM  1 , Incompressible 
X Planar  Tension,  Po,  N=2,  RM  1 , Incompressible 
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Figure  4.15  - Calibration  to  planar  tension  without  compressibility  using  RM1 
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Engineering  Tensile  Stress  (MPa)  Engineering  Tensile  Stress  (MPa) 


Planar  Stress  found  from  VL  function 

O Planar  Tension,  Og,  N=6,  VL1 , Compressible 
A Planar  Tension,  Og,  N=3,  VL1,  Compressible 
x Planar  Tension,  Po,  N=2,  VL1,  Compressible 


16  - Calibration  to  planar  tension  with  compressibility  using  VL1 


Planar  Stress  found  from  VL  function 

o Planar  Tension,  Og,  N=6,  VL2,  Compressible 
A Planar  Tension,  Og,  N=3,  VL2,  Compressible 
x Planar  Tension,  Po,  N=2,  VL2,  Compressible 


Figure  4.17  - Calibration  to  planar  tension  with  compressibility  using  VL2 
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□ Based  on  Volumetric  Compression  Test 
O Volumetric  Compression,  Og,  N=6,  RM 1 
X Volumetric  Compression,  Og,  N=6,  VL1 

160 
140 
120 
100 
80 
60 
40 
20 
0 

0.95  0.96  0.97  0.98  0.99  1 

Volume  Ratio 

Figure  4.18  - Calibration  to  volumetric  compression  test  using  RM1  and  VL1  models 


Reduced  Tensile  Stress  (MPa).  Experimental 

♦—  Reduced  Compressive  Stress  (MPa),  Experimental 
o Reduced  Uniaxial  Stress,  Og6,  based  on  RM1 

a Reduced  Uniaxial  Stress,  Og3,  based  on  RM  1 

x Reduced  Uniaxial  Stress,  Pol2,  based  on  RM1 
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Figure  4.19  - Reduced  uniaxial  stresses  from  experiment  vs  ABAQUS  results,  using 
RM1 
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Reduced  Tensile  Stress  (MPa),  Experimental 

♦—  Reduced  Compressive  Stress  (MPa),  Experimental 
o Reduced  Uniaxial  Stress,  Og6,  based  on  VL1 

a Reduced  Uniaxial  Stress,  Og3,  based  on  VL1 

x Reduced  Uniaxial  Stress,  Pol2,  based  on  VL1 
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Figure  4.20  - Reduced  uniaxial  stresses  from  experiment  vs  ABAQUS  results,  using 
VL1 
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Reduced  Tensile  Stress  (MPa),  Experimental 

♦—  Reduced  Compressive  Stress  (MPa),  Experimental 
o Reduced  Uniaxial  Stress,  Og6,  based  on  VL2 

a Reduced  Uniaxial  Stress,  Og3,  based  on  VL2 

x Reduced  Uniaxial  Stress,  Pol2,  based  on  VL2 
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Figure  4.21  - Reduced  uniaxial  stresses  from  experiment  vs.  ABAQUS  results,  using 
VL2 
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5.  FINITE  ELEMENT  MODELING  OF  STEEL 


When  compressively  loading  the  bearings,  the  steel  shims  undergo  plastic 
straining,  and  eventually  some  steel  shims  rupture,  causing  the  bearing  to  lose  load- 
carrying  capacity.  Therefore,  in  order  to  predict  the  shape  of  the  bearing  compressive 
stress-strain  curve  into  the  failure  region,  the  plastic  deformation  of  the  steel  must  be 
modeled.  Uniaxial  tensile  tests  on  steel  specimens  have  shown  a linear  elastic 
response  followed  by  a strain-hardening  plastic  range.  This  type  of  response  can  be 
modeled  in  most  finite  element  codes.  In  ABAQUS,  the  material’s  elastic  behavior  is 
defined  by  furnishing  the  Elastic  Modulus  and  Poisson’s  ratio  (0.3  in  all  runs)  under 
the  ELASTIC  option.  The  plastic  data  are  defined  in  the  PLASTIC  option  by  giving 
the  current  yield  stress  as  a function  of  the  plastic  strain. 

Prior  to  shim  failure,  necking  occurs.  Necking  is  localized  yielding  in  the 
steel.  A triaxial  state  of  stress  exists  during  necking.  The  model  used  in  this  analysis 
cannot  accurately  predict  the  plastic  flow  at  this  point.  Consequently,  the  steel  model 
cannot  predict  failure  of  the  shims.  However,  we  are  able  to  capture  the  gross  plastic 
flow  and  therefore  get  an  accurate  depiction  of  the  deformed  shape  of  the  bearing  for 
much  of  the  loading  history. 

5.1  Steel  Testing  Methods 

For  the  B-Bearings,  the  shims  were  made  of  ASTM  A570  Grade  36  steel.  To 
prepare  the  shims  for  the  bonding  process  with  the  rubber,  they  are  sandblasted  and 
then  coated  with  a chemical  compound  in  a process  known  as  chemlocking.  Tests 
were  done  on  smooth  coupons,  sandblasted  coupons,  and  chemlocked  coupons  in 
order  to  determine  if  sandblasting  and  chemlocking  had  a noticeable  effect  on  the  steel 
material  properties.  A change  in  behavior  by  sandblasting  and/or  chemlocking  must 
be  taken  into  account  so  that  accurate  predictions  of  full-scale  bearing  behavior  can  be 
made  from  scale  model  behavior.  For  instance,  in  sandblasting,  fine  grains  of  sand  are 
propelled  at  the  specimen  under  high  pressure,  giving  the  surface  a pitted,  uneven 
texture.  The  depth  of  this  unevenness  would  be  a higher  percentage  of  the  overall 
thickness  for  a thinner  steel  specimen,  and  may  change  its  mechanical  properties. 
Chemlock,  on  the  other  hand,  is  an  adhesive  that  is  painted  on  the  sandblasted 
specimen  in  order  to  improve  the  steel/rubber  bond  during  the  vulcanization  process. 
Since  the  chemlock  is  applied  only  on  the  surface  of  the  steel,  it  appears  unlikely  that 
it  would  have  any  significant  affect  on  the  steel  properties. 

Five  steel  sheets  were  provided  from  which  dumbbell  specimens  were  cut, 
three  specimens  from  each  of  the  five  steel  sheets,  for  a total  of  fifteen  tensile 
specimens,  as  shown  in  Table  5.1. 
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Table  5.1  - Tensile  dogbone  specimens  by  finish 


Sheet  Thickness 

Finish 

# tensile  specimens 

Specimen 

Designation 

12-gage 

smooth 

3 

S-#12  * 

12-gage 

sandblasted 

3 

SB-#12 

12-gage 

chemlocked 

3 

CL-#12 

16-gage 

sandblasted 

3 

SB-#16 

22-gage 

sandblasted 

3 

SB-#22 

* # refers  to  specimen  1,  2,  or  3. 

The  tension  tests  were  conducted  in  accordance  with  ASTM  E 8-95a,  Standard 
Test  Methods  for  Tension  Testing  of  Metallic  Materials.  Figure  5.1  shows  the 
dimensions  for  the  tensile  specimens.  The  tolerance  on  the  width  of  the  reduced 
section  is  per  the  ASTM  Standards,  as  is  the  minimum  radius  on  the  fillets.  The  ends 
of  the  specimen  were  cut  at  least  3.00  inches  in  length  to  ensure  adequate  gripping  in 
our  test  machine. 

The  relevant  stress  and  strain  measures  are  true  (Cauchy)  stress  and 
logarithmic  strain  for  input  into  ABAQUS,  and  material  data  for  all  these  models  must 
be  given  in  these  measures.  The  testing  machine  automatically  generated  an 
engineering  stress-strain  curve  for  each  test.  Special  steps  were  necessary  to  find  the 
true  stress  and  logarithmic  strain,  and  are  described  below. 

During  plastic  straining  and  prior  to  necking,  it  is  well  established  that  a metal 
maintains  its  volume,  i.e.  it  is  essentially  incompressible  (Dowling,  1993;  Dieter, 
1976).  This  observation  allows  a simple  transformation  from  engineering  stress  (er)  to 

A l 

true  stress  (/).  Since  the  material  is  incompressible,  Al  = A0l0  =>  — = — , where  l 

A l0 

is  the  current  gage  length,  /0is  the  original  gage  length,  and  Aa  is  the  original  cross- 
sectional  area.  Therefore, 


_ F 
“ A 


F_ 

4 


= <7- 


— = cr(l  + e) 

h 


where  e is  engineering  strain. 


Also,  the  logarithmic  strain  £ is  defined  as 


(5.1) 
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e - In 


= ln(l  + e) 


(5.2) 


Strictly  speaking.  Equation  (5.1)  only  applies  in  the  plastic  range  up  to 
necking.  This  equation  was  applied  to  the  elastic  region  also,  since  the  error  in  doing 
so  is  no  larger  than  e , which  is  0.002  at  yielding.  The  correct  true  stress  falls  between 
the  engineering  stress  and  the  true  stress  as  given  by  Equation  (5.1),  so  0.2%  is  an 
upper  bound  on  the  error. 

Once  necking  starts,  nearly  all  of  the  strain  is  localized  at  the  neck.  Since 
Equations  (5.1)  and  (5.2)  represent  an  average  strain  across  the  gage  length,  they  are 
not  valid.  Actual  measurement  of  the  cross-sectional  area  at  the  point  of  necking  must 


be  made,  and  t - — and  e = In 
A 1 


The  width  of  the  neck  was  measured  using  a caliper  with  a digital  readout 
accurate  to  0.01  mm.  Three  micrometers  were  used  to  measure  the  thickness  of  the 
specimens  prior  to  loading.  Figure  5.2  shows  the  micrometer  designs. 

Micrometer  #1,  with  circular  anvils,  had  a vernier  which  allowed 
measurements  to  an  accuracy  of  0.0001  in.  Micrometers  #2  and  #3  were  accurate  to 
0.001  in.  Micrometer  #1  was  the  most  accurate  instrument,  and  the  greatest 
confidence  was  placed  in  its  measurements.  Each  micrometer  was  used  to  measure 
the  original  thickness  in  order  to  ensure  that  measurements  between  them  correlated 
well,  which  they  did.  However,  the  circular  anvils  of  Micrometer  #1  were  too  large  to 
measure  the  thickness  of  the  specimens  in  the  necking  region.  This  was  also  true  for 
Micrometer  #2.  Micrometer  #3’s  anvils  came  to  a point,  allowing  more  accurate 
measurements  in  the  necking  region.  Therefore,  Micrometer  #3  was  used  to  measure 
the  thickness  in  the  necking  region. 

The  thickness  measurement  in  the  necking  region  was  tricky.  First,  at  the 
onset  of  necking,  which  occurs  at  the  maximum  force  F (Dowling,  1993;  Dieter, 

1976),  necking  was  not  visible.  Once  the  load  started  to  drop,  the  test  machine  was 
stopped  at  intervals,  and  the  necking  measurements  were  made.  It  was  noticed  that  the 
location  where  minimum  measurements  were  recorded  shortly  after  the  maximum 
force  F was  reached  (in  which  necking  was  not  visible)  did  in  fact  correspond  to  the 
location  where  necking  became  visible  shortly  thereafter.  Therefore,  once  localized 
straining  occurred,  it  didn’t  shift  to  another  location,  so  the  location  of  the  initial 
minimum  measurements  corresponded  with  the  location  of  the  later  minimum 
measurements.  Second,  the  neck  wasn’t  perpendicular  to  the  longitudinal  axis.  If  the 
failure  mechanism  is  shearing,  or  sliding  of  the  molecules,  the  maximum  shear  should 

be  at  a 45°  angle  to  the  longitudinal  axis,  which  is  the  axis  of  loading.  In  practice,  this 
angle  varied  between  90°  and  45° , sometimes  fairly  close  to  90° , but  in  general  closer 
to  45° . Care  had  to  be  taken  to  make  the  measurements  along  this  line. 
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5.2  Steel  Testing  Results 


Figures  5.3  through  5.7  show  the  results  of  each  of  the  five  sets  of  tests. 
Necking  begins  when  each  of  the  solid  lines  ends.  The  stress-strain  points  in  this 
region  are  based  on  actual  measurements  of  thickness  in  the  necking  region.  The 
shape  of  the  curve  in  the  necking  region  is  generally  linear,  which  is  in  close 
agreement  with  true  stress-strain  curves  shown  in  Dowling  (1993)  and  Dieter  (1976). 
There  is  some  variation,  and  this  is  primarily  due  to  inaccurate  recording  of  the  load  at 
which  the  test  machine  was  stopped  (in  order  to  take  the  thickness  measurements)  and 
inaccurate  thickness  measurements.  These  variations  from  linearity  are  more 
prevalent  at  higher  strains,  due  in  some  cases  to  the  fact  that  the  load  was  rapidly 
declining  at  higher  strains,  leading  to  inaccurate  recording  of  load. 

The  repeatability  of  each  set  of  tests  was  good.  The  stresses  for  each  of  the 
three  curves  for  each  sample  differ  by  less  than  1.5%  in  the  plastic  range,  with  the 
exception  of  the  22-gage  tests,  where  the  stresses  for  SB3-22  are  about  7%  larger  than 
SB  1-22.  There  are  only  two  specimens  in  the  SB#-12  series  (one  specimen  in  the 
series  was  damaged  during  testing). 

Figure  5.8  is  a comparison  by  finish  of  the  three  12-gage  samples,  which 
demonstrates  that  the  stress  for  the  chemlocked  specimens  was  significantly  lower  in 
the  plastic  region  than  for  the  smooth  or  sandblasted  specimens.  We  requested  the 
bearing  manufacturer  to  supply  us  with  steel  samples  from  the  same  millings  used  in 
the  manufacture  of  the  shims.  However,  we  cannot  be  sure  that  the  samples  received 
were  from  the  same  millings.  Differences  in  steel  behavior  may  be  due  to  the  steel 
coming  from  different  steel  millings.  It  is  common  to  see  yield  stresses  between  248 
MPa  (36  ksi)  and  345  MPa  (50  ksi)  from  milling  to  milling  for  the  steel  used  in  the  B- 
Bearings,  ASTM  A570  Grade  36  steel.  Nevertheless,  in  our  samples,  the  yield 
stresses  between  the  12-gage  samples  were  260  MPa  for  the  sandblasted  and 
chemlocked  specimens  versus  270  MPa  for  the  smooth  specimens,  a difference  of 
only  about  3.8%.  The  larger  difference  in  behavior  shows  up  in  the  plastic  range, 
where  the  stresses,  for  strains  larger  than  0.2,  are  about  13%  lower  for  the  chemlocked 
specimen  compared  to  the  sandblasted  and  smooth  specimens. 

Figure  5.9  is  a comparison  between  the  12, 16,  and  22-gage  specimens,  and 
shows  that  the  stress  for  CL1-12,  over  the  plastic  range  of  0.2  to  1.0,  is  within  1.0%  of 
SB1-16  and  no  more  than  6.0%  less  than  SB2-22. 

5.3  Measurement  Accuracy  and  Variability  and  Effect  on  Bearing  Response 

A Satec  tension  and  compression  machine  was  used  to  test  the  steel  specimens. 
The  accuracy  of  this  machine  is  established  through  a calibration  using  known  dead 
loads,  and  was  calibrated  according  to  ASTM  E4,  which  requires  that  the  machine 

reading  be  within  ± 1%  of  the  true  reading.  The  machine  has  met  these  requirements. 
This  is  equivalent  to  a relative  expanded  uncertainty  (k=2)  of  1%,  and  is  small 
compared  to  specimen  variability. 

Table  5.2  summarizes  the  testing  results  for  the  elastic  properties.  Note  that 
the  elastic  moduli  are  about  the  same  for  each  type  of  specimen,  but  the  yield  stress  is 
about  10%  higher  for  the  22-gage  specimen. 
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Table  5.2  - Comparison  of  elastic  properties  for  steel  specimens 


Steel  Specimen 

Yield  Stress 
(MPa) 

Elastic  Modulus 
(E) 

12-gage  Smooth  Finish  (S3- 12) 

270 

2.1x10s 

12-gage  Sandblasted  Finish  (SB2-12) 

262 

2.1x10s 

12-gage  Chemlocked  Finish  (CL1-12) 

262 

2.2  x10s 

16-gage  Sandblasted  Finish  (SB  1-16) 

261 

2.3  x10s 

22-gage  Sandblasted  Finish  (SB2-22) 

290 

2.0  x10s 

The  test  results  for  each  specimen  in  the  plastic  range  fall  into  two  distinct 
groups  in  which  the  stress-strain  responses  are  similar.  In  the  first  group,  consisting  of 
the  12-gage  smooth  and  sandblasted  specimens,  the  stresses  are  more  than  10%  higher 
than  the  stresses  in  the  second  group,  consisting  of  the  12-gage  chemlocked  specimen, 
the  16-gage  chemlocked  specimen,  and  the  22-gage  sandblasted  specimen  (c.f.  Figures 
5.8  and  5.9).  Figure  5. 10  compares  the  second  group  more  closely,  up  to  a logarithmic 
strain  of  0.60,  and  shows  that  the  stresses  of  the  22-gage  sandblasted  shim  are 
approximately  4%  to  6%  higher  than  the  stresses  of  the  12-gage  chemlocked  specimen 
and  the  16-gage  sandblasted  specimen.  The  differences  in  these  samples  indicate 
differences  in  material  properties,  since  the  variations  are  larger  than  the  resolution  of 
the  testing  machine. 

Based  on  these  testing  results,  we  expect  similar  bearing  behavior  predicted  by 
the  finite  element  analysis  up  to  initial  shim  yielding  when  incorporating  any  of  the 
12-gage  samples  or  the  16-gage  sample  in  the  model,  but  higher  stresses  at  initial  shim 
yielding  when  incorporating  the  22-gage  sample  (assuming  usage  of  the  same  rubber 
model  and  model  geometry  for  all  runs).  From  initial  shim  yielding  to  compressive 
strains  near  bearing  failure,  the  steel  is  yielding,  and  we  expect  the  predicted  response 
from  the  finite  element  analysis  to  be  virtually  the  same  with  inclusion  of  either  the 
12-gage  chemlocked  sampleor  the  16-gage  sandblasted  sample,  and  approximately  4% 
to  6%  higher  for  the  22-gage  sandblasted  sample  (c.f.  Figures  5.9  and  5.10).  We 
expect  the  stresses  in  the  plastic  range  to  be  significantly  higher  than  the  foregoing 
with  inclusion  of  either  the  12-gage  smooth  or  sandblasted  specimens  in  the  model 
(c.f.  Figure  5.8).  Finite  element  analysis  of  the  bearings,  discussed  in  detail  in  Section 
6,  verify  the  validity  of  the  expected  responses. 
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Figure  5.1  - Dimensions  of  tensile  dogbone  specimens 
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Figure  5.2  - Micrometers  used  to  measure  the  necking  region  of  the  tensile  dogbone 
specimens 
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5.3  - True  stress-strain  curves  for  smooth  (S)  12-gage  specimens 
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Figure  5.4  - True  stress-strain  curves  for  sandblasted  (SB)  12-gage  specimens 
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Figure  5.5  - True  stress-strain  curves  for  chemlocked  (CL)  12-gage  specimens 
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Figure  5.6  - True  stress-strain  curves  for  sandblasted  (SB)  16-gage  specimens 
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Figure  5.7  - True  stress-strain  Curves  for  sandblasted  (SB)  22-gage  specimens  by 
finish 
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Figure  5.8  - Comparison  of  true  stress-strain  curves  for  12-gage  steel  samples 
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Figure  5.9  - Comparison  of  12,  16,  and  22-gage  steel  samples 
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Figure  5.10  - Comparison  of  12,  16,  and  22-gage  steel  samples  up  to  a logarithmic 
strain  of  0.60 


88 


6.  FINITE  ELEMENT  MODELING  OF  THE  ELASTOMERIC  BEARINGS 

6.1  Discretization  of  the  Finite  Element  Models  of  the  Bearings 

As  noted  in  Section  5,  the  finite  element  model  does  not  have  the  capability  to 
model  failure  of  the  steel  shims,  since  the  plastic  flow  of  the  steel  in  the  necking 
region  cannot  be  accurately  modeled.  This  is  not  a drawback  in  the  context  of  this 
research,  because  we  are  primarily  interested  in  predicting  the  ultimate  capacity  of  a 
full-scale  bearing  based  on  a test  to  failure  of  a scale  model  bearing.  Comparison  of 
stress-strain  curves  between  scale  models  and  prototypes  will  be  the  primary  means  of 
determining  whether  the  failure  stress  of  a scale  model  is  representative  of  a full-scale 
bearing.  The  results  of  scale  model  analyses  are  discussed  in  Section  6.5. 

Assuming  the  load  to  be  uniformly  distributed  over  the  bearings,  an 
axisymmetric  model  was  used.  There  is  also  symmetry  between  the  top  and  bottom 
half  of  the  bearing.  Therefore,  only  1/4  of  the  bearing  was  modeled.  The  model 
represents  a slice  of  the  bearing  through  the  R-Z  plane,  from  R = 0 to  the  outer  edge, 
and  from  the  middle  to  the  top. 

Each  rubber  layer  was  modeled  with  four  layers  of  2-dimensional  linear 
displacement  continuum  axisymmetric  hybrid  elements  with  full  integration 
(CAX4H),  and  each  steel  shim  was  modeled  with  CAX4  elements,  the  same  as  for 
rubber  except  the  hybrid  feature  is  not  used.  Since  the  rubber  at  the  cover  makes 
contact  with  the  endplate  in  the  course  of  the  loading  history,  the  endplate  and  cover 
were  modeled  using  the  ABAQUS  CONTACT  PAIR  option  with  the  endplate  as  the 
master  surface  and  the  outside  surfaces  of  certain  elements  in  the  cover  as  the  slave 
surface.  Since  deformations  of  the  endplate  are  negligible,  it  was  modeled  using  the 
1 -dimensional  rigid  elements  RAX2.  The  static  coefficient  of  friction  between  the 
rubber  and  end  plate  is  ji  = 0.5,  based  on  a simple  test  which  determined  the  critical 

friction  angle  as  45°.  Figures  6.1  through  6.3  show  the  undeformed  meshes  for  the 
Full-Scale,  1/2-Scale,  and  1/4-Scale  models,  respectively.  All  deformed  mesh  and 
contour  plots  were  generated  in  PATRAN.  The  X-Y  plane  in  these  PATRAN 
generated  plots  corresponds  to  the  R-Z  plane  of  the  axisymmetric  model. 

6.2  Effect  of  Different  Factors  on  the  Results  of  the  Finite  Element  Model 
6.2.1  Inclusion  of  the  Central  Core  in  the  Model 

After  testing  to  failure,  the  B-Bearings  were  cut  open.  Figure  6.4  is  a 
photograph  of  a vertical  slice  through  the  core  of  Bearing  B2,  the  1/2-Scale  Bearing. 
This  photograph  provides  a baseline  on  which  to  judge  the  accuracy  of  the  deformed 
shapes  generated  by  the  finite  element  analysis. 

As  noted  earlier,  the  core  of  the  bearing  is  filled  with  rubber  after  stacking  of 
the  rubber  and  shim  layers.  We  attempted  to  model  the  bearing  with  the  core.  The 
ABAQUS  filename  for  this  model  is  procore_RMl_CS7_CLl  12:1,  which  means 
Prototype  with  the  core  using  the  RM1  rubber  model,  volumetric  data  CS7,  and  steel 
data  CL1-12.  However,  the  analysis  failed  to  converge  beyond  a bearing  compressive 
strain  of  14.3%,  which  is  54%  of  the  strain  when  the  actual  bearing  failed.  Figure  6.5 
shows  the  deformation  in  the  rubber  and  steel  at  this  point.  We  believe  the  analytical 
problem  is  due  to  excessive  distortion  of  the  rubber  elements  in  the  core  that  are 
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adjacent  to  the  steel  elements  (all  elements  were  initially  rectangular  in  shape).  In  the 
actual  bearing,  the  stresses  in  the  core  rubber  at  the  interface  with  the  bonded  portion 
of  the  bearing  are  so  great  that  the  core  rubber  splits  away  from  the  rubber  and  steel  in 
this  region.  This  would  have  the  affect  of  reducing  the  stresses  in  this  region.  This 
loss  of  contact  has  not  been  modeled. 

By  excluding  the  core,  we  were  able  to  get  the  analyses  to  go  out  beyond  the 
bearing  compressive  strains  at  which  bearing  failure  occured  (finite  element  run 
entitled  “pro_RMl_CS7_CL  112:1”).  Figure  6.6  shows  the  ABAQUS  runs  with  and 
without  the  core  in  comparison  with  the  experimental  Prototype  Bearing  stress-strain 
curve.  Based  on  this  result,  it  is  reasonable  to  assume  that  the  rubber  core  has 
negligible  effect  on  the  behavior  of  the  bearing.  Furthermore,  the  finite  element  model 
performs  well  in  estimating  the  experimental  curve,  with  the  exception  that  the  finite 
element  model  is  stiffer  prior  to  yielding  of  the  shims.  The  abrupt  change  in  modulus 
of  the  bearing  stress-strain  curve  is  due  to  yielding  of  the  shims.  For  the  experiment 
on  the  Prototype,  this  occured  at  approximately  42  MPa,  versus  about  40  MPa  for  the 
finite  element  run  pro_RMl_CS7_CLl  12:1.  A secant  modulus  at  the  point  of  initial 
shim  yielding  gives  667  MPa  for  the  experimental  test  versus  1025  MPa  for  the  finite 
element  run  pro_RMl_CS7_CLl  12:1.  A straight  line  drawn  from  the  shim  yielding 
point  to  failure  gives  235  MPa  for  the  experimental  test  versus  215  MPa  for  the  finite 
element  run  pro_RMl_CS7_CLl  12: 1 . Failure  for  the  bearing  occurs  when  a shim 
ruptures,  which  decreases  the  ability  of  the  bearing  to  sustain  load,  and  appears  on  the 
bearing  stress-strain  curve  as  a sudden  drop  at  a bearing  engineering  compressive 
stress  of  90  MPa  and  strain  of  0.267  m/m.  The  finite  element  run 
pro_RMl_CS7_CLl  12:1  gives  a stress  of  88.9  MPa  at  a strain  of  0.267  m/m,  or  1.2% 
lower  than  experiment.  Finite  element  runs  with  CL1-12  will  be  used  as  the  baseline 
for  comparisons  between  models. 

6.2.2  Comparison  between  the  Ogden  and  Polynomial  Strain  Energy  Functions 
and  Rubber  Models  RM1  and  VL1 

The  finite  element  runs  plotted  in  Figure  6.6,  procore_RMl_CS7_CLl  12:1 
and  pro_RMl_CS7_CLl  12:1,  were  conducted  using  the  Ogden  Strain  Energy 
Potential  with  N=6.  Figure  6.7  is  a comparison  of  pro_RMl_CS7_CLl  12: 1 with  a 
finite  element  run  using  the  Polynomial  Strain  Energy  Function,  N=2  (filename 
pro_RMl_CS7_CLl  12_Po2:l).  The  results  are  virtually  identical.  All  other  runs 
discussed  in  the  following  sections  use  the  Ogden  Strain  Energy  Potential  with  N=6. 

Figure  6.8  shows  the  engineering  stress-strain  curve  for  the  Prototype  Bearing, 
along  with  the  finite  element  results  using  rubber  models  RM1  and  VL1,  the 
volumetric  data  from  CompSet  7 (CS7),  and  steel  data  CL  1-1 2.  Both  models  give  the 
same  response.  We  are  tempted  to  attribute  this  to  the  fact  that  the  finite  element 
program,  using  RM1  (which  does  not  include  the  planar  tensile  stresses),  accurately 
estimates  the  planar  stresses  found  from  the  Valanis-Landel  function,  and  that  are 
included  in  VL1.  In  other  words,  RM1  and  VL1  are  virtually  identical.  However, 
these  rubber  models  only  contribute  to  the  deviatoric  portion  of  the  strain  energy 
function.  To  reinforce  this  statement,  the  Ogden  Strain  Energy  Potential  is  rewritten 
here. 
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The  first  summation  of  the  strain  energy  potential  is  the  deviatoric  strain 
energy.  Rubber  models  RM1  and  VL1  are  used  to  find  only  the  parameters  and  ctj, 
which  define  the  deviatoric  strain  energy.  The  parameter  Z),  is  found  from  the 
volumetric  test  data  (CS7),  and  determines  the  volumetric  strain  energy.  Both  models 
use  CS7  to  determine  the  volumetric  strain  energy. 

As  shown  in  Section  6.3,  the  hydrostatic  stresses  are  two  orders  of  magnitude 
greater  than  the  deviatoric  stresses,  so  that  if  RM1  and  VL1  had  yielded  different 
deviatoric  stresses,  the  compressive  response  of  the  bearing  would  have  been  affected 
very  little. 

6.2.3  Rubber  Compressibility 

The  next  two  figures.  Figures  6.9  and  6.10,  look  at  the  effect  of  compressibility 
on  the  results.  An  incompressible  rubber  model  makes  the  bearing  much  stiffer  prior 
to  yielding  of  the  shims.  The  initial  compressive  modulus  for  the  bearing  modeled 
with  the  incompressible  rubber  assumption  (pro_RMl_Incomp_CLl  12:1)  is  1745 
MPa,  which  is  70%  larger  (1025  MPa)  than  die  initial  compressive  modulus  of  the 
analysis  done  which  included  compressibility  of  the  rubber 
(pro_RMl_CS7_CLl  12:1).  Beyond  initial  shim  yielding,  although  the  predicted 
stresses  are  higher  in  the  incompressible  model,  the  stiffness,  based  on  the  slope  of  a 
straight  line  drawn  from  the  point  of  initial  shim  yielding  to  the  strain  at  which  the 
bearing  failed,  is  virtually  the  same  (214  MPa  for  the  Incompressible  model  versus 
215  MPa  for  the  Compressible  model).  For  the  incompressible  model,  the  stress  at  the 
strain  at  which  the  bearing  failed  is  90.3  MPa,  which  is  only  0.3%  higher  than  the 
stress  at  which  the  actual  bearing  failed  (90  MPa),  and  only  1%  higher  than  the  failure 
stress  in  the  compressible  model  (88.9  MPa).  Therefore,  the  inclusion  of 
compressibility  in  the  rubber  model  has  a large  effect  on  the  compressive  stiffness  of 
the  bearing  prior  to  initial  shim  yielding,  but  has  a small  effect  on  the  predicted  stress 
at  failure. 

6.2.4  Steel  Test  Data 

Figure  6. 1 1 compares  the  finite  element  model  using  the  chemlocked  steel  data 
versus  the  finite  element  model  using  the  sandblasted  steel  data.  The  response  for  the 
bearing  model  based  on  steel  specimen  SB2-12  is  stiffer  in  the  plastic  range  compared 
to  steel  specimen  CL1-12,  and  this  is  due  to  the  difference  in  the  tensile  specimen 
stress-strain  response  of  these  steel  models  in  the  plastic  range,  as  shown  in  Figure 
6.12.  Using  specimen  SB2-12,  the  bearing  engineering  compressive  stress  at  the 
strain  at  which  the  bearing  failed  (0.267  m/m)  is  99.7  MPa,  or  10.8%  higher  than 
experiment.  Although  not  apparent  from  Figure  6.26,  due  to  the  large  scale  shown  for 
the  logarithmic  strain,  the  yield  stresses  for  the  two  steel  models  are  nearly  identical 
(261.9  MPa  for  steel  specimen  CL1-12  versus  261.6  MPa  for  steel  specimen  SB2-12). 
This  is  reflected  in  the  identical  shape  of  the  bearing  stress-strain  curves  up  to  initial 
yielding  of  the  shims,  as  shown  in  Figure  6.1 1. 

6.3  Deformation  and  Contour  Plots  and  the  Dominance  of  Hydrostatic  Stress 

Figure  6.13  is  a plot  of  the  finite  element  solution  of  the  deformed  shape  of  the 
rubber  portion  of  the  bearing  for  model  pro_RMl_CS7_CLl  12:1  at  Increment  96, 
which  corresponds  to  a bearing  engineering  stress  of  88.6  MPa  and  a bearing 
engineering  strain  of  0.265  m/m.  Consulting  Figure  6.8,  this  is  very  near  the 
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compressive  engineering  strain  at  which  the  Prototype  failed.  The  deformed  shape  of 
the  rubber  is  very  similar  to  the  actual  deformed  shape  of  the  1/2-Scale  Bearing 
(Figure  6.4).  Even  though  we  are  showing  only  the  rubber  elements  here,  the 
deformed  shape  of  the  shims  is  clearly  visible  between  the  rubber  layers,  and  is  also 
close  to  the  actual  deformed  shape  of  the  shims  of  the  1/2-Scale  Bearing. 

Superimposed  on  the  bearing  is  the  hydrostatic  stress  in  the  rubber,  which 
reaches  a maximum  of  83.7  MPa.  The  hydrostatic  stress  is  an  invariant,  and  is  the 
average  of  the  three  normal  stresses.  Figures  6.14  through  6.16  show  the  three 
principal  stresses,  each  of  which  has  a maximum  value  of  about  83  MPa.  Since  the 
principal  stresses  have  approximately  the  same  magnitude,  the  maximum  shearing 
stresses  must  be  relatively  small.  Therefore,  we  conclude  that  for  the  loading 
condition  modeled,  the  hydrostatic  stresses  are  much  larger  than  the  deviatoric 
stresses. 

The  only  plane  in  which  shearing  can  occur  is  the  R-Z  plane,  due  to  symmetry. 
Figure  6.17  shows  the  maximum  shearing  stresses.  For  presentation  of  scalar  data, 
which  is  shown  in  contour  plots,  PATRAN  takes  the  largest  and  smallest  values  of  the 
scalar  variable,  and  divides  the  spectrum  evenly.  In  Figure  6.17,  the  largest  value  of 
the  maximum  shear  is  65.7  MPa,  but  this  occurs  only  in  a very  few  elements  near  the 
cover  layer  of  the  bearing.  The  great  majority  of  the  elements  have  a maximum  shear 
stress  between  zero  and  7.30  MPa.  To  better  see  the  distribution  of  stresses  in  this 
lowest  range  of  the  spectrum,  Figure  6. 1 8 is  shown,  where  the  largest  value  was 
limited  to  7.30  MPa.  This  figure  clearly  shows  that  the  majority  of  the  elements  have 
a maximum  shearing  stress  less  than  1 .46  MPa.  Of  course,  any  element  in  the  highest 
spectrum  range  has  a maximum  shear  value  at  least  as  large  as  that  range,  but  this  is  a 
relatively  small  set  of  elements  which  are  concentrated  near  the  cover  layer.  The 
relatively  low  maximum  shear  values  are  not  surprising,  since  rubber  is  weak  in  shear. 

Figure  6.19  is  the  finite  element  solution  of  the  deformed  shape  of  the  rubber 
portion  of  the  bearing,  along  with  the  hydrostatic  stresses,  for  the  same  model 
(pro_RMl_CS7_CLl  12: 1)  at  Increment  8,  which  corresponds  to  a bearing 
engineering  stress  of  39.5  MPa  and  a bearing  engineering  strain  of  0.0383  m/m.  This 
is  the  location  on  the  engineering  stress-strain  curve  where  the  slope  changes  abruptly, 
indicating  yielding  of  the  shims.  At  Increment  8,  the  hydrostatic  stresses  fall  off  more 
rapidly  as  we  move  toward  the  outer  edge  of  the  bearing  as  compared  to  Increment  96. 
The  maximum  value  of  the  hydrostatic  stress  is  68.7  MPa,  or  about  82%  of  the 
maximum  hydrostatic  stress  at  Increment  96. 

Figures  6.20  through  6.22  show  the  three  principal  stresses  in  the  rubber  at 
Increment  8,  and  Figure  6.23  shows  the  maximum  shearing  stress  in  this  same 
increment.  Similar  to  Increment  96,  the  hydrostatic  stress  is  much  larger  than  the 
maximum  shearing  stress;  e.g.,  from  the  core  to  the  middle  of  the  bearing,  the 
hydrostatic  stress  is  between  53.0  and  68.7  MPa,  and  the  maximum  shearing  stress  is 
between  zero  and  0.792  MPa. 

Figures  6.24  and  6.25  show  the  von  Mises  stresses  in  the  steel  at  Increments  8 
and  96,  respectively.  As  noted  above,  when  modeling  the  core  of  the  bearing,  we  were 
not  able  to  get  convergence  beyond  14.3%  compressive  engineering  strain,  most  likely 
due  to  the  severe  distortion  of  the  rubber  elements  adjacent  to  the  core.  Therefore,  we 
neglected  the  core  and  prevented  movement  in  the  R-direction  of  all  elements,  rubber 
and  steel,  that  are  adjacent  to  the  core.  Therefore,  the  stresses  and  strains  near  the  core 


92 


will  overstate  the  actual  stresses  and  strains  in  the  bearing.  In  the  actual  bearing,  the 
shims  are  not  constrained  against  movement  in  the  R-direction,  and  the  core  rubber 
separates  from  the  shims,  relieving  stresses  on  the  shims  along  the  core.  Shim  failure 
in  the  actual  bearing  does  not  occur  at  the  core,  but  actually  occurs  less  than  halfway 
from  the  core  to  the  cover  layer.  Figure  6.26  shows  the  equivalent  plastic  strain  at 
Increment  96.  The  highest  plastic  strains  are  at  the  interior  of  the  shim  near  the 
midheight  of  the  bearing,  as  we  would  expect.  In  spite  of  these  inaccurate  boundary 
conditions,  as  we  move  farther  from  the  core  the  stresses  and  strains  should  more 
closely  approximate  the  actual  values.  This  is  clear  when  we  look  at  the  deformed 
shape  of  the  bearing.  The  simulated  shim  deformation  accurately  predicts  the  actual 
shim  deformation  as  we  move  toward  the  cover  layer  of  the  bearing. 

6.4  Discussion  of  Differences  in  Initial  Stiffness  of  Finite  Element  Model  versus 
Experiment 

Figure  6.27  shows  the  stress-strain  curve  in  the  region  prior  to  the  abrupt 
change  in  the  stiffness.  Why  is  the  finite  element  model  significantly  stiffer  than 
experiment?  In  the  experimental  results,  low  stiffness  at  low  stresses  is  present  in 
each  of  the  bearings.  This  initial  low  stiffness  is  due  in  part  or  entirely  to  plate  and 
head  misalignment,  and/or  the  endplates  being  out  of  parallel.  This  effect  has  not  been 
modeled.  Other  possible  causes  of  the  initial  low  stiffness  may  include  prestresses  of 
the  rubber  and  steel  due  to  the  vulcanization  process,  and  prestresses  of  the  steel  due 
to  cold  rolling. 

6.5  Results  of  Finite  Element  Analysis  on  Scale  Models  and  Geometric 
Dissimilarities 

The  actual  geometry  of  the  1/2  and  1/4-Scale  models  was  reproduced  in  the 
finite  element  models,  including  the  geometric  dissimilarities  discussed  in  Section 
2.1.1.  Figure  6.28  shows  the  1/2-Scale  experimental  bearing  stress-strain  curve,  and 
the  finite  element  models  using  the  compressible  rubber  model  (volumetric  data  CS7) 
with  two  different  steel  models,  the  first  model  based  on  the  16-gage  sandblasted  steel 
specimen  SB  1-1 6,  and  the  second  model  based  on  the  12-gage  chemlocked  steel 
specimen  CL1-12.  The  responses  of  the  finite  element  models  are  virtually  the  same, 
reflecting  the  similarity  of  die  steel  stress-strain  curves,  as  shown  in  Figure  6.29.  The 
yield  stresses  for  the  steel  models  are  virtually  identical  (261.9  MPa  for  steel 
specimen  CL  1-1 2 versus  260.5  MPa  for  specimen  SB  1-16).  This  is  reflected  in  the 
nearly  identical  shape  of  the  bearing  stress-strain  curve  up  to  initial  yielding  of  the 
shims,  as  shown  in  Figure  6.28.  The  actual  bearing  failed  at  a compressive  strain  of 
0.227  m/m  and  a compressive  stress  of  88.7  MPa.  The  model  based  on  steel  specimen 
SB  1-1 6 (l/2Scale_RMl_CS7_SBl  16)  has  a compressive  stress  of  90.3  MPa  at  the 
failure  strain  of  0.227,  and  for  the  model  based  on  steel  specimen  CL1-12  the 
compressive  stress  at  the  failure  strain  is  90.8  MPa  (l/2-Scale_RMl_CS7_CLl  12) 

The  secant  modulus  at  the  point  of  initial  shim  yielding  is  730  MPa  for  the  actual  test 
versus  1033  MPa  for  both  finite  element  models.  The  models  overestimated  the 
bearing  stiffness  in  this  region  by  41.5%.  The  slope  of  a straight  line  drawn  from 
initial  shim  yielding  to  bearing  failure  for  the  actual  test  is  263  MPa  versus  252  MPa 
for  l/2-Scale_RMl_CS7_CLl  12  and  249  MPa  for  l/2-Scale_RMl_CS7_SBl  16. 

Figure  6.30  shows  the  1/4-Scale  experimental  bearing  stress-strain  curve,  and 
the  finite  element  models  using  the  compressible  rubber  model  (volumetric  data  CS7) 
with  two  different  steel  models,  the  first  model  based  on  the  22-gage  sandblasted  steel 
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specimen  SB2-22,  and  the  second  model  based  on  the  12-gage  chemlocked  steel 
specimen  CL1-12.  The  response  of  the  finite  element  models  reflects  the  differences 
of  the  steel  stress-strain  curves,  as  shown  in  Figure  6.31.  The  yield  stress  for  steel 
specimen  SB2-22  is  289.9  MPa,  or  10.7%  larger  than  the  yield  stress  of  261.9  MPa  for 
specimen  CL1 12.  The  shims  yielded  in  the  bearing  test  at  approximately  51.5  MPa, 
but  for  l/4-Scale_RMl_CS7_CLl  12  initial  yielding  occured  at  42.9  MPa  and  for  1/4- 
Scale_RMl_CS7_SB222  at  55.3  MPa.  The  actual  bearing  failed  at  a compressive 
strain  of  0.28 1 m/m  and  a compressive  stress  of  96.7  MPa.  Model 
1/4  Scale_RMl_CS7_SB222  has  a compressive  stress  of  105  MPa  at  the  failure  strain 
of  0.281,  and  in  model  l/4-Scale_RMl_CS7_CLl  12  the  compressive  stress  at  the 
failure  strain  is  100.4  MPa.  The  secant  modulus  at  the  point  of  initial  shim  yielding  is 
472  MPa  for  the  actual  test  versus  approximately  975  MPa  for  both  finite  element 
models.  The  finite  element  models  overestimated  the  bearing  stiffness  in  this  region 
by  207%.  However,  as  noted  in  Section  2.6.6,  the  tangent  stiffness  of  the 
experimental  curve  is  930  MPa  just  prior  to  the  initial  yielding  of  the  shims,  which  is 
very  close  to  the  slope  of  each  finite  element  run.  This  is  apparent  in  Figure  6.30. 

The  initial  low  stiffness  at  low  stresses  is  present  in  each  of  the  experimental  bearings 
tested,  but  is  particularly  severe  for  the  1/4-Scale  Bearing.  As  discussed  in  Section 
6.3,  this  initial  low  stiffness  may  be  due  to  plate  and  head  misalignment,  and/or  the 
end  plates  being  out  of  parallel.  These  affects  would  be  magnified  for  smaller  scale 
bearings,  since  they  are  made  to  the  same  tolerances  as  full-scale  bearings.  The  slope 
of  a straight  line  drawn  from  initial  shim  yielding  to  bearing  failure  for  the  actual  test 
is  263  MPa  versus  223  MPa  for  1/4  Scale_RMl_CS7_SB222  and  242  MPa  for 
1/4  Scale_RMl_CS7_CLl  12.  The  finite  element  results  are  accurate,  as  is  apparent 
by  inspection  of  Figure  6.30,  especially  after  taking  into  account  the  initial  low 
stiffness  of  the  actual  1 /4-Scale  Bearing. 

As  noted  in  Section  2.1.1,  the  B-Bearings  are  dissimilar  in  three  ways:  shim 
thickness,  center  alignment  hole  diameter,  and  cover  layer  thickness.  The  1/4  and  1/2- 
Scale  shims  are  22-gage  and  16-gage,  respectively,  or  nominally  0.762  mm  (0.03  in) 
and  1 .524  mm  (0.06  in)  thick,  and  are  therefore  similar  in  thickness,  but  the  Prototype 
shim  is  12-gage,  or  nominally  2.667  mm  thick  (0.105  in),  so  it  is  0.381  mm  (0.015  in) 
thinner,  or  12.5%  thinner,  than  it  should  be  in  order  to  have  a geometrically  similar 
thickness.  The  diameter  of  the  center  alignment  hole  for  the  Prototype  is  50.1  mm,  but 
the  diameters  of  the  center  alignment  holes  for  the  1/2-Scale  and  1/4-Scale  models  are 
each  30.1  mm.  The  1/2-Scale  model's  center  alignment  hole  diameter  is  20.1%  larger 
than  it  should  have  been  (25.05  mm)  for  geometric  similarity,  and  the  1 /4-Scale 
model's  center  alignment  hole  diameter  is  141%  larger  than  it  should  have  been 
(12.5025  mm)  for  geometric  similarity.  The  cover  layer  for  the  Prototype  is  19  mm 
thick,  and  should  be  9.5  mm  thick  for  the  1/2-Scale  model,  but  is  actually  12.5  mm,  or 
32%  too  large,  and  should  be  4.75  mm  thick  for  the  1/4-Scale  model,  but  is  actually  19 
mm,  or  300%  too  thick  for  strict  geometric  similarity  to  hold.  For  convenience,  this 
information  is  compiled  in  Table  6.1  shown  immediately  below. 
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Table  6. 1 - Summary  of  geometric  dissimilarities 


Prototype  (mm) 

1/2  Scale  (mm) 

1/4  Scale  (mm) 

Shim  Thickness 

Actual  Dimensions 

2.667 

1.524 

0.762 

For  Geometric  Similarity 

3.048 

1.524 

0.762 

Center  Alignment  Hole 

Diameter 

Actual  Dimensions 

50.1 

30.1 

30.1 

For  Geometric  Similarity 

50.1 

25.05 

12.5025 

Cover  Layer  Thickness 

Actual  Dimensions 

19.0 

12.5 

19.0 

For  Geometric  Similarity 

19.0 

9.50 

4.75 

A finite  element  run  was  conducted  for  the  Prototype  using  steel  model  CL1  - 
12,  deviatoric  rubber  model  RM1,  and  volumetric  test  data  CS7,  but  with  the  shims 
made  geometrically  similar  to  the  1/2  and  1/4-Scale  models;  i.e.  3.048  mm  thick  (11 
gage  shims).  Figure  6.32  shows  the  comparison  of  this  run  with  each  of  the  bearings 
as-is.  Using  geometrically  similar  shims  for  the  prototype,  the  figure  shows  a large 
increase  in  the  stress-strain  response  of  the  bearing  between  initial  shim  yielding  and 
failure.  For  example,  the  stress  at  a strain  of  0.25  m/m  was  10%  higher  using  1 1 gage 
shims.  Furthermore,  the  response  of  the  prototype,  using  similar  shims,  is  nearly 
identical  to  the  1/4  and  1/2-Scale  Bearing  responses. 

The  geometric  dissimilarities  in  the  central  alignment  hole  and  cover  layer  are 
most  pronounced  with  the  1/4-Scale  Bearing,  as  shown  in  Table  6.1.  As  shown  in 
Figure  6.32,  in  spite  of  these  geometric  dissimilarities,  the  response  of  the  1 /4-Scale 
Bearing  is  nearly  identical  to  the  prototype,  assuming  geometrically  similar  shims. 
Figure  6.33  shows  the  separate  effects  of  each  of  the  dissimilarities,  and  they  are  seen, 
in  this  particular  case,  to  cancel  each  other  out.  For  example,  making  the  core  of  the 
1 /4-Scale  Bearing  geometrically  similar  to  the  prototype  (leaving  all  other  geometry 
unchanged)  increases  the  stress  by  approximately  3%  at  a strain  of  0.25  m/m,  and 
making  the  cover  layer  of  the  1 /4-Scale  Bearing  geometrically  similar  to  the  prototype 
(leaving  all  other  geometry  unchanged)  decreases  the  stress  by  about  the  same  amount 
at  the  same  strain  of  0.25  m/m.  Clearly,  the  effects  of  the  dissimilarities  of  the  core 
and  cover  layer  on  the  response  of  the  1/4-Scale  Bearing  are  less  important  than  the 
dissimilarity  in  the  shim  thickness. 
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Figure  6.1-  Undeformed  mesh  for  Full-Scale  Bearing 


96 


Figure  6.2  - Undeformed  mesh  for  1/2-Scale  Bearing 
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Figure  6.3  - Undeformed  mesh  for  1/4-Scale  Bearing 
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Figure  6.4  - Photograph  of  central  slice  of  failed  1/2-Scale  Bearing  (B2) 
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Figure  6.5  - Deformed  shape  and  steel  plastic  strains  in  the  prototype  modeled  with 

the  central  rubber  core 
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Figure  6.6  - Comparison  for  the  prototype  between  the  experimental  engineering  stress 
strain  curve  and  ABAQUS  finite  element  runs  with  and  without  the 
modeling  of  the  core 
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Figure  6.7  - Comparison  of  ABAQUS  runs  based  on  Ogden,  N=6  and  Polynomial, 
N=2  Strain  Energy  Functions 
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Figure  6.8  - Comparison  for  the  prototype  between  the  experimental  engineering  stress 
strain  curve  and  the  ABAQUS  finite  element  runs  based  on  RM1  and 
VL1 
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Figure  6.9  - Comparison  for  the  prototype  between  the  experimental  engineering 
stress-strain  curve  and  the  ABAQUS  finite  element  runs  using  a 
compressible  versus  an  incompressible  rubber  model 
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Figure  6. 10  - Low  strain  region  comparison  for  the  prototype  between  the 

experimental  engineering  stress-strain  curve  and  the  ABAQUS  finite 
element  runs  using  a compressible  versus  an  incompressible  rubber 
model 
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Figure  6.1 1 - Comparison  for  the  prototype  between  the  experimental  engineering 
stress-strain  curve  and  the  ABAQUS  finite  element  runs  using  the 
chemlocked  steel  data  (CL1-12)  versus  the  sandblasted  steel  data  (SB2- 
12) 
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Figure  6.12  - Comparison  between  the  true  stress-strain  curves  of  the  chemlocked  steel 
data  (CL  1-12)  versus  the  sandblasted  steel  data  (SB2-12) 
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MSC/PATRAN  Version  6.0  06-Aug-97  15:23:12 

FRINGE:  Static,  Step1,TotalTime=  1.4405:  Stress,  Components  (HYDROSTATIC)  -ABAQUS 
DEFORMATION:  Static,  Step1,TotalTime=1.4405:  Deformation,  Displacements  -ABAQUS 


Figure  6.13  - Full-Scale  Bearing  - Deformation  and  hydrostatic  stresses  in  the  rubber 
at  the  compressive  engineering  failure  strain 
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Figure  6. 14  - Full-Scale  Bearing  - Deformation  and  major  principal  stresses  in  the 
rubber  at  the  compressive  engineering  failure  strain 
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Figure  6.15  - Full-Scale  Bearing  - Deformation  and  intermediate  principal  stresses  in 
the  rubber  at  the  compressive  engineering  failure  strain 
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Figure  6.16  - Full-Scale  Bearing  - Deformation  and  minor  principal  stresses  in  the 
rubber  at  the  compressive  engineering  failure  strain 
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Figure  6.17  - Full-Scale  Bearing  - Deformation  and  maximum  shearing  stresses  in  the 
rubber  at  the  compressive  engineering  failure  strain 
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Figure  6.18  - Full-Scale  Bearing  - Deformation  and  maximum  shearing  stresses  (up  to 
7.30  MPa)  in  the  rubber  at  the  compressive  engineering  failure  strain 
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Figure  6.19  - Full-Scale  Bearing  - Deformation  and  hydrostatic  stresses  in  the  rubber 
at  the  compressive  engineering  strain  corresponding  to  initial  shim 
yielding 
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Figure  6.20  - Full-Scale  Bearing  - Deformation  and  major  principal  stresses  in  the 
rubber  at  the  compressive  engineering  strain  corresponding  to  initial 
shim  yielding 
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Figure  6.21  - Full-Scale  Bearing  - Deformation  and  intermediate  principal  stresses  in 
the  rubber  at  the  compressive  engineering  strain  corresponding  to  initial 
shim  yielding 
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Figure  6.22  - Full-Scale  Bearing  - Deformation  and  minor  principal  stresses  in  the 
rubber  at  the  compressive  engineering  strain  corresponding  to  initial 
shim  yielding 
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Figure  6.23  - Full-Scale  Bearing  - Deformation  and  maximum  shearing  stresses  in  the 
rubber  at  the  compressive  engineering  strain  corresponding  to  initial 
shim  yielding 
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Figure  6.24  - Full-Scale  Bearing  - Deformation  and  von  Mises  stresses  in  the  steel  at 
the  compressive  engineering  strain  corresponding  to  initial  shim 
yielding 
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Figure  6.25  - Full-Scale  Bearing  - Deformation  and  von  Mises  stresses  in  the  steel  at 
the  compressive  engineering  failure  strain 
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Figure  6.26  - Full-Scale  Bearing  - Deformation  and  equivalent  plastic  strain  in  the 
steel  at  the  compressive  failure  engineering  strain 
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Figure  6.27  - Low  strain  region  comparison  for  the  prototype  between  the 

experimental  engineering  stress-strain  curve  and  the  ABAQUS  finite 
element  runs  based  on  RM 1 and  VL 
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Figure  6.28  - Comparison  for  the  1 /2-Scale  Bearing  model  between  the  experimental 
engineering  stress-strain  curve  and  the  ABAQUS  finite  element  runs 
using  the  chemlocked  steel  data  (CL  1-1 2)  versus  the  sandblasted  steel 
data  (SB  1-16) 
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Figure  6.29  - Comparison  between  the  true  stress-strain  curves  of  the  chemlocked  steel 
data  (CL  1-12)  versus  the  sandblasted  steel  data  (SB  1-16) 
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Figure  6.30  - Comparison  for  the  1 /4-Scale  Bearing  model  between  the  experimental 
engineering  stress-strain  curve  and  the  ABAQUS  finite  element  runs 
using  the  chemlocked  steel  data  (CL1-12)  versus  the  sandblasted  steel 
data  (SB2-22) 
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Figure  6.31  - Comparison  between  the  true  stress-strain  curves  of  the  chemiocked  steel 
data  (CL1-12)  versus  the  sandblasted  steel  data  (SB2-22) 
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Figure  6.32  - Comparison  between  the  ABAQUS  finite  element  runs  for  the  prototype, 
1/2  and  1/4  scale  models  using  geometrically  nonsimilar  and  similar 
shims  for  the  prototype 
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Figure  6.33  - Separate  and  combined  effects  of  core  and  cover  layer  dissimilarities  on 
the  response  of  the  1/4-Scale  Bearing 
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7.  CONCLUSIONS 


Control  of  structural  vibrations  through  the  use  of  elastomeric  bearings  is  a 
very  promising  technology.  It  is  important  to  have  a good  understanding  of  the  loads 
that  may  cause  failure  of  elastomeric  bearings,  so  that  a design  engineer  may  have  an 
accurate  assessment  of  the  factor  of  safety  built  into  his  or  her  design.  As  a step  in 
this  direction,  in  this  study  we  have  documented  failure  loads  of  elastomeric  bearings 
in  compression,  in  order  to  determine  if  scale  models  can  be  used  to  predict  the  failure 
load  of  full-scale  bearings.  Testing  of  full-scale  elastomeric  bearings  to  failure  in 
compression  is  not  possible  in  most  testing  facilities  due  to  the  limitation  of  loading 
capacity.  Even  though  the  elastomeric  bearings  were  designed  to  be  exactly 
geometrically  similar,  certain  dissimilarities  were  present  due  to  manufacturing  errors 
and  to  the  use  of  standard  manufacturing  bearing  molds.  These  dissimilarities  have 
been  analyzed  using  finite  element  models  to  determine  their  influence  on  the  overall 
performance  of  the  bearings. 

The  findings  and  conclusions  of  this  study  are  summarized  below: 

1 . In  order  to  predict  failure  loads  of  full  scale  bearings  from  scale  models,  the  stress- 
strain  responses  of  full-scale  and  scale  models  must  be  nearly  identical.  Assuming 
identical  rubber  and  steel  properties  among  all  bearings,  the  use  of  geometrically 
similar  shims  and  rubber  layers  resulted  in  virtually  identical  stress-strain 
responses  for  all  bearings,  in  spite  of  the  dissimilarities  in  the  central  alignment 
hole  and  cover  layer.  The  use  of  a 12  gage  shim  (2.667  mm  thick)  instead  of  a 
geometrically  similar  1 1 gage  shim  (3.048  mm  thick)  for  the  Full-Scale  Bearing 
reduced  the  stresses  significantly  (10%  reduction  at  a strain  of  0.25  m/m). 
Therefore,  geometric  similarity  in  shim  thickness  and  rubber  thickness  is  most 
important.  In  actual  practice,  failure  stresses  may  vary  due  to  material 
imperfections  and  imperfect  manufacture. 

2.  The  effects  of  the  dissimilarities  in  the  scale  models  in  the  central  alignment  holes 
and  the  cover  layers  are  less  important  than  shim  and  rubber  thickness,  and  in  this 
particular  case  cancel  each  other  out.  The  dissimilarities  are  most  prevalent  for  the 
1 /4-Scale  Bearing,  and  using  a similar  central  alignment  hole  (everything  else 
unchanged)  for  the  1 /4-Scale  Bearing  increases  the  stress  by  approximately  3%  at 
a strain  of  0.25  m/m,  and  using  a similar  cover  layer  (everything  else  unchanged) 
decreases  the  stress  by  the  same  amount  at  0.25  m/m.  For  the  Full-Scale  Bearing, 
the  diameter  of  the  center  alignment  hole  was  50. 1 mm,  or  6.6%  of  the  bonded 
diameter.  The  diameter  of  the  central  alignment  hole  for  the  1 /4-Scale  Bearing 
was  141%  larger  than  it  should  have  been  in  order  to  be  geometrically  similar  to 
the  Full-Scale  Bearing,  resulting  in  a smaller  bonded  area  than  if  it  were 
geometrically  similar.  The  cover  layer  for  the  Full-Scale  Bearing  was  19  mm 
thick,  or  5.0%  of  the  bonded  radius.  The  cover  layer  for  the  1 /4-Scale  Bearing  was 
300%  larger  than  it  should  have  been  in  order  to  be  geometrically  similar  to  the 
Full-Scale  Bearing. 

3.  In  the  experimental  bearing  tests,  the  1/2-Scale  Bearing  compressive  stress-strain 
curve  is  stiffer  than  the  Full-Scale's  curve  in  the  plastic  range  of  the  bearing,  i.e. 
between  initial  shim  yielding  and  bearing  failure.  This  difference  in  stiffness  can 
be  attributed  to  the  fact  that  the  shims  of  the  Full-Scale  Bearing  are  less  than  twice 
the  thickness  of  the  1/2-Scale  Bearing.  The  finite  element  analysis  shows  that 
bearings  with  geometrically  similar  shims  and  identical  steel  properties  have 
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identical  compressive  stress-strain  responses.  However,  the  experimental  1/4- 
Scale  Bearing  stress-strain  response  is  stiffer  than  the  1 /2-Scale  Bearing  response 
in  the  bearing's  plastic  region,  and  shows  a higher  stress  at  initial  shim  yielding,  in 
spite  of  similarity  in  shim  thickness  between  these  models.  Tensile  tests  on  22- 
gage  sandblasted  coupons,  which  were  used  in  the  1/4-Scale  Bearing,  showed  a 
yield  stress  approximately  10%  higher  than  the  16-gage  sandblasted  samples  (used 
in  the  1/2-Scale  Bearing),  and  stresses  in  the  plastic  range  approximately  4%  to 
6%  higher.  The  shims  for  our  1/4-Scale  Bearing  (22-gage)  are  only  0.762  mm 
thick,  much  thinner  than  normally  used  in  bearing  construction.  The  cold  rolling 
process  used  in  the  manufacturing  of  this  sheet  steel  and  the  scouring  of  the  metal 
surface  due  to  sandblasting  may  have  modified  its  mechanical  properties. 

4.  The  Valanis-Landel  function  was  used  to  extend  the  experimental  rubber  stress- 
strain  data  (uniaxial  tension  and  compression)  to  the  planar  tensile  (pure  shear) 
mode.  The  planar  tensile  stress  values  are  slightly  higher  than  the  uniaxial  tensile 
stresses,  which  agrees  with  experimental  data  reported  by  Treloar  (1943). 
Furthermore,  equibiaxial  stresses  found  from  the  conversion  of  uniaxial 
compressive  data  agree  with  the  equibiaxial  stresses  found  from  the  Valanis- 
Landel  function.  Using  only  uniaxial  data  to  determine  the  elastic  constants,  the 
Ogden  Strain  Energy  Potential,  which  is  a special  case  of  the  Valanis-Landel 
function,  is  able  to  predict  the  planar  tensile  stresses,  while  the  Polynomial  Strain 
Energy  Potential  with  second  order  terms  (N=2)  overestimates  the  planar  tensile 
stresses.  While  the  qualitative  agreement  with  Treloar’ s experimental  data  is 
compelling,  the  definitive  proof  of  the  ability  of  the  Valanis-Landel  function  to 
accurately  predict  planar  stresses  from  uniaxial  stresses  requires  a comparison  with 
experimental  planar  stress  data.  In  the  future,  it  would  be  interesting  to  compare 
experimental  planar  stresses  to  those  predicted  by  the  Valanis-Landel  function 
from  only  uniaxial  data. 

5.  The  finite  element  results  on  the  bearings  suggest  that  the  compressive  stress-strain 
response  of  the  bearings,  for  a given  steel  model,  is  dependent  primarily  on  the 
bulk  modulus  of  rubber,  and  for  finite  element  studies  where  the  hydrostatic 
rubber  stress  dominates,  experimental  determination  of  the  bulk  modulus  may  be 
all  that  is  necessary  to  adequately  represent  the  strain  energy  function.  Stress- 
strain  data  in  the  uniaxial,  equibiaxial,  and  planar  deformation  modes  were  used  to 
find  the  deviatoric  parameters  of  the  strain  energy  function.  The  hydrostatic 
parameters  were  found  from  the  volumetric  test  data  and  were  used  throughout  the 
analysis  when  compressibility  was  incorporated  in  the  model,  hence  no  differences 
were  observed  in  the  hydrostatic  pressure. 

6.  In  the  actual  bearings,  the  central  alignment  hole  is  filled  with  rubber.  Including 
this  central  core  of  rubber  in  the  finite  element  model  resulted  in  convergence 
problems,  most  likely  due  to  excessive  deformation  of  the  rubber  elements  in  the 
core  adjacent  to  the  shims.  Excluding  the  core  resulted  in  no  measurable 
difference  in  the  stress-strain  response  of  the  bearing.  The  core  represents  0.4%  of 
the  bonded  area  (including  the  core)  for  the  Full-Scale  Bearing,  0.6%  for  the  1/2- 
Scale  Bearing,  and  2.5%  for  the  1/4-Scale  Bearing. 

7.  The  inclusion  of  a compressible  rubber  model  had  little  affect  on  the  predicted 
failure  stress,  even  though  the  initial  compressive  modulus  was  much  closer  to  the 
compressive  modulus  in  the  experiments.  The  stress  at  the  strain  at  which  the 
Full-Scale  Bearing  failed  was  90.3  MPa  for  the  finite  element  model  based  on 
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rubber  incompressibility,  which  is  only  1%  larger  than  the  failure  stress  with  the 
inclusion  of  compressibility  in  the  finite  element  model  (88.9  MPa).  The  Full- 
Scale  Bearing  failed  at  90  MPa,  which  falls  between  the  incompressible  and 
compressible  finite  element  results.  The  initial  compressive  modulus  of  the  model 
of  the  Full-Scale  Bearing  was  1745  MPa  when  an  incompressible  rubber  model 
was  used,  which  is  70%  larger  than  the  initial  compressive  modulus  when  a 
compressible  rubber  model  was  used  (1025  MPa).  The  shape  factor  of  the  Full- 
Scale  Bearing  rubber  layers  was  15  (same  for  scale  models),  and  the  thickness  of 
the  12  gage  shims  (2.67mm)  was  21%  of  the  thickness  of  the  rubber  layers  (12.7 
mm  or  1/2  in.). 
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APPENDIX  A 


SOME  ELEMENTARY  RESULTS  IN  LARGE  STRAIN  KINEMATICS 

The  following  brief  description  of  large  strain  kinematics  is  based  largely  on 
Hjelmstad’s  “Fundamentals  of  Structural  Mechanics”  (1997).  The  detailed 
implementation  of  the  kinematics  and  constitutive  theory  used  for  rubber  and  steel 
modeling  in  ABAQUS  can  be  found  in  the  ABAQUS/Standard  Users  Manual  (1995a) 
and  the  ABAQUS  Theory  Manual  (1995b). 

In  one  dimension,  the  stretch  A = j-  is  the  ratio  of  the  deformed  length  of  a 

line  to  the  original  length  of  the  line.  In  terms  of  the  stretch,  familiar  definitions  of 
strain  are 
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More  generally,  the  stretch  can  vary  from  point  to  point  in  a one-dimensional 

^0 

rod,  so  that  a rod  of  initial  length  4>,  when  deformed,  has  a length  / = J A(x)dx.  If  the 

o 

strain  in  the  rod  is  homogeneous,  i.e.  the  same  at  each  point,  then 

'o 

/ = A J dx  - 


A = 


/ 


To  get  a general  formulation  for  the  stretch  at  a point  in  three  dimensions, 

consider  Figure  A.l.  Curve  C is  any  curve  in  the  undeformed  body  that  passes 

through  point  P,  and  is  deformed  into  the  line  O(C).  If  we  parametrize  C by  s,  the 

X(s  + As)-  X(s)^\  dX  dX.  D,.  . n D 

— ' v ' -=,  and is  a unit  vector  at  P that  is  tangent  to  C at  P, 

ds  ds 


lim 

As-»0' 


As 


7 


because  the  length  of  the  vector  X(s  + As)-  X(s)  approaches  the  arclength  s as  As 


approaches  zero.  The  jimj 


x(s  + A sy 
As  y 


dx, 

= — is  a vector  tangent  to  O(C)  at  O(P),  but 
ds 
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is  not  a unit  vector.  The  stretch  at  point  P in  the  direction  given  by  parameter  s along 

~ . - dx 

curve C is  A-  — =. 

ds 

Undeformed  Body  Deformed  Body 


Figure  A.l  - Mapping  from  the  Undeformed  to  Deformed  Configuration 

There  are  an  infinite  number  of  lines  that  pass  through  the  point  P,  each 
defining  a direction  with  an  associated  stretch,  so  the  state  of  deformation  at  P can  be 
thought  of  as  a infinite  collection  of  stretches  and  their  associated  directions.  To  get  a 
general  description  which  captures  all  of  the  information,  we  can  relate  the  vector 
dx/ds  to  dX/ds.  Writing  x(s)  in  terms  of  the  deformed  body’s  coordinate  system,  we 
get  x(j)  = jc,(5)  i,  + jc20)  i2  + Jt3(.y)  /3 . In  terms  of  the  mapping  function  F,  this  can  be 
expressed  as 


x(s)  = 0,(X(s))  h + 02(X(j))  i2  + <J>3(X(s))  i3 
Taking  the  derivative  of  Equation  (A1.6)  with  respect  to  s,  we  get 


(A  1.6) 
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Dropping  the  reference  to  the  arclength  parameter  s,  and  putting  Equation 
(A  1.7)  in  matrix  form,  we  get 
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F is  called  the  “deformation  gradient,”  and  is  a tensor  of  rank  two.  By  the 
polar  decomposition  theorem,  F=RU,  where  U changes  the  length  of  a vector 
without  changing  its  direction,  and  R rotates  the  vector  without  changing  its  length. 
The  Right  Cauchy-Green  Deformation  Tensor  is  defined  as  C = FT  F,  which  is  by 
definition  symmetric.  C contains  information  on  deformation  only,  as  can  be  shown 

by  C = FrF  = (Ruf  RU  = UTRTRU  = UTU.  Since  R is  a rotation  only,  / = RT  R. 


The  Lagrangian  Strain  in  three  dimensions  is 


e=\(£-D- 


Since  C is  symmetric,  E 


is  also  symmetric.  In  indicial  notation,  Ei}  = ui}  + uj{  +uk  iuk  j. 


In  terms  of  principal  stretches, 


^ 0 0 > 

0/^0 


The  invariants  of  C are 


/,  = AJ  + A^  + A^ 

I2  = AjA^  + A^A^  + A^ 


h = WA 


131 


The  ratio  of  the  current  volume  to  the  original  volume  is 
J - = det(F)  = AjA^Aj . As  noted  in  Section  3.2,  in  the  mixed  formulation  revised 

invariants  are  used  to  remove  the  volume  change.  In  ABAQUS,  the  revised 
invariants,  written  in  terms  of  the  deviatoric  stretches.  A,  = J~V3Xi , are 

72  = ^ A^  + Af  A^  + A^A| 

As  shown  in  Section  3.2,  the  constitutive  models  for  rubber  are  written  in 
terms  of  the  strain  invariants  or  principal  stretches.  ABAQUS  provides  output  in 
terms  of  strain.  For  steel,  the  strain  is  decomposed  into  an  elastic  and  inelastic 
portion. 
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APPENDIX  B 


PROGRAM  LISTING  FOR  DETERMINATION  OF  THE  DERIVATIVE  OF 
THE  VALANIS-LANDEL  FUNCTION 

/*  This  program  uses  the  Kearsley  and  Zappas  separation  algorithm 
to  calculate  the  Valanis-Landel  function  from  uniaxial  compression 
and  tension  data. 

*/ 

#include  <stdio.h> 

#include  <math.h> 

double  SeriesTerm(double*); 

int  main(void){ 

enum  {Max=100}; 


int  i=0,j=l,iCount=0; 

double  dArgl=0.,  dArg2=0.,  dLamda=0.,  dExpl=0.,  dExp2=0.; 
double  dTotal,  dTotal_old,  dF,  dF_old,  dG,  dG_old; 
double  dVL[2]; 
double  dTol  = le-04; 
double  dTruel,  dTrue2; 

double  *dpArgl,  *dpArg2; 

printf("\n\nTolerance  = %f\n\n",  dTol); 

printf("%-4s  %-5s  %-\  Is  %-10s  %-10s  %-10s\n\n", 

"Max", "#Iters", "Stretch", "Sum[i]",  "Sum[i-1] "/'Difference"); 

for(dLamda  = 0.15;  dLamda  <=  6.5;  dLamda  +=  0.20)  { 

dTotal  = dTotal_old  = 0.0; 
dF  = dF_old  = 0.0; 
dG  = dG_old  = 0.0; 
i=l; 

for(i=l;i<=Max;i++)  { 
dExpl  = pow(0.25,i-l); 
dExp2  = -0.5  * pow(0.25,i-l); 
dArgl  = pow(dLamda, dExpl); 
dpArgl  = &dArgl; 
dArg2  = pow(dLamda,dExp2); 
dpArg2  = &dArg2; 
dF_old  = dF; 
dG_old  = dG; 
dTotal_old  = dF  + dG; 

dF  = dF  + SeriesTerm(dpArgl); 
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dG  = dG  + SeriesTerm(dpArg2); 
dTotal  = dF  + dG; 

dVL[0]  = dTotal_old/dLamda; 
dVL[l]  = dTotal/dLamda; 
if(fabs(dVL[0]  - dVL[l])  < dTol)  break; 

} /*  End  of  Iteration  Loop  */ 

printf("%-7d  %-5d  %-8.2f  %-l  1.4f  %-12.4f  %-10.4f\n".  Max,  i-1, 
dLamda,dVL[l],dVL[0],  dVL[l]  - dVL[0]); 

} /*  End  of  Loop  that  finds  VL  function  values  (at  Lamda)  */ 

} /*  End  of  program  */ 


double  SeriesTerm(double  * dpArg)  { 

double  dTerm; 
double  dArg  = *dpArg; 


if(dArg<=0.3)  { 

dTerm  = -1.1 6648 17e02  + 9.8062739e02  * dArg 

+ -4.5649689e02  * pow(dArg,2)  + -2.1791899e04  * pow(dArg,3) 
+ 8.4299544e04  * pow(dArg,4)  + -9.6543082e04  * pow(dArg,5); 

return  dTerm; 

} 

else  if(dArg>0.3  &&  dArg<=0.55)  { 

dTerm  = -3.4861 124e01  + 2.4064177e02  * dArg 

+ -6.9647254e02  * pow(dArg,2)  + 9.763 1903e02  * pow(dArg,3) 
+ -6.1623463e02  * pow(dArg,4)  + 1.1560896e02  * pow(dArg,5); 

return  dTerm; 


else  if(dArg>0.55  &&  dArg<=0.70)  { 

dTerm  = 1.246665e01  + -1. 390855 le02  * dArg 

+ 4.9764865e02  * pow(dArg,2)  + -8.4097151e02  * pow(dArg,3) 
+ 6.9239043e02  * pow(dArg,4)  + -2.2408201e02  * pow(dArg,5); 

return  dTerm; 

} 

else  if(dArg>0.70  &&  dArg<=0.80)  { 
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dTerm  = -8.8530545  + 2.7500424e01  * dArg 

+ -3.334608 led  * pow(dArg,2)  + 1.4650992e01  * pow(dArg,3); 

return  dTerm; 

} 

else  if(dArg>0.80  &&  dArg<=0.90)  { 

dTerm  = 1.4608343e03  + -8.741 1359e03  * dArg 

+ 2.0866959e04  * pow(dArg,2)  + -2.48636  le04  * pow(dArg,3) 

+ 1. 478989 le04  * pow(dArg,4)  + -3.5133156e03  * pow(dArg,5); 

return  dTerm; 

} 

else  if(dArg>0.90  &&  dArg<=0.95)  { 

dTerm  = 3.247632 le03  + -1.4074399e04  * dArg 

+ 2.2859003e04  * pow(dArg,2)  + -1.6496004e04  * pow(dArg,3) 
+ 4.463750 le03  * pow(dArg,4); 

return  dTerm; 

} 

else  if(dArg>0.95  &&  dArg<=0.997 1 1 ) { 

dTerm  = 1 .274908  le03  + -4.8341207e03  * dArg 

+ 5.513 1936e03  * pow(dArg,2)  + -1.4714537e02  * pow(dArg,3) 
+ -3.295223e03  * pow(dArg,4)  + 1.488383e03  * pow(dArg,5); 

return  dTerm; 

} 

else  if(dArg>0.997 1 1 &&  dArg<=1.0)  { 
dTerm  = 12.216955  * (dArg  - 1); 
return  dTerm; 

} 

else  if(dArg>1.0  &&  dArg<=  1.0274)  { 
dTerm  = 14.254380  * (dArg  - 1); 
return  dTerm; 

} 

else  if(dArg> 1.0274  &&  dArg<=  1.0621)  { 
dTerm  = 5.03458  * dArg  - 4.78196; 
return  dTerm; 
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else  if(dArg> 1.0621  &&  dArg<=2.4971)  { 

dTerm  = -1.9199173e03  + 9.8937267e03  * dArg 

+ -2.235642  le04  * pow(dArg,2)  + 2.9064259e04  * pow(dArg,3) 

+ -2.395601 9e04  * pow(dArg,4)  + 1.2989462e04  * pow(dArg,5) 

+ -4.6372447e03  * pow(dArg,6)  + 1.0523308e03  * pow(dArg,7) 

+ -1.3795277e02  * pow(dArg,8)  + 7.9740265  * pow(dArg,9); 
return  dTerm; 

} 

else  { 

dTerm  = -1.8322249e03  + 4.8664372e03  * dArg 

+ -5.5502795e03  * pow(dArg,2)  + 3.5819692e03  * pow(dArg,3) 

+ -1. 444874  le03  * pow(dArg,4)  + 3.7870 137e02  * pow(dArg,5) 

+ -6.4597472e01  * pow(dArg,6)  + 6.9249393  * pow(dArg,7) 

+ -4.2390496e-01  * pow(dArg,8)  + 1.1303178e-02  * pow(dArg,9); 
return  dTerm; 

} 

} 
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